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STATICAL CHARGES OF DEFORMED POLYMER 


By CzEstaw Bazan 
Physical Laboratory of B. Bierut University, Wroclaw 
(Received May 15, 1958) 


Thin plates of a polyvinylchloride clamped from two sides, heated above the softening 
point and then cooled to room temperature acquire electrical charges. These charges have 
opposite sign on both sides of the sample. The plates prepared in this way show optical 
double refraction. From the observed patterns one can conclude that the chain segments 
are oriented preferentially in the radial direction. This effect may be related to the electri- 
zation of the sample. It is suggested that the assymetry of the electrical charges may be 
due to polar group orientation. 


§ 1. Introduction 


It is known that a deformed polymer exhibits optical double refraction. This 
effect is caused by the orientation of the chain segments, which are aligned paralelly. 
Such an orientation may occur under pressing, stretching or by viscous flow in solution. 
At yet this effect has not been related to static electrization of the polymer materials. 
Skinner (1953) explains the formation of such static charges by charge carriers migra- 
ting into the polymer from the surrounding material (especially metal). This mechanism 
is not connected with orientation effects in the polymer. In this paper it is shown 
that a correlation between electrization and orientation exists. 


§ 2. Experimental results i 
a. Mechanical properties 


An investigation was made of a pure polyvinychloride. Discs were formed from 

the polymer under a pressure of 330 atm. They were white, opaque and showed no 

` electrical charge. Then they were heated over 3 hours at a temperature of 90—100" G 
that is, above the softening point, and afterwards cooled to room temperature, During 

the heating process the sample was clamped between aluminium, tin or copper foil, 

or between bakelite plates. The pressure was of the order 15 atm. In this way one 

obtains transparent samples, 45 mm in diameter and 1,1 mm thick, with a density 


of 1,36 g.cm™. f 
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After reheating the sample in the same manner, but without pressure, the diameter 
of the sample diminishes and the thickness increases; these changes are of the order 
10%. The density is constant within 1%. 

After repeated heating under clamping the deformation occurs. The deformation 
process is thus reversible, in agreement with observations of Kargin and Sokolova 
(1957). In time the sample takes on a yellow hue, which indicates some irreversible 
changes. 

The observed changes of the dimensions of the sample demonstrate that a rever- 
sible “flow of molecules in the radial direction takes place in the polymer. This flow 
is caused by the thermal motion of the chain segments above the softening point 
(Buchdahl 1950, Floyd 1952). The applied pressure tends to align these segments in 
radial direction. When the sample is heated without pressure, the termal agitation 
destroys this alignement. 


b. Optical properties 


Double refraction in the deformed and undeformed polymer was examined 
by observing them in linearly polarised light between crossed polarizers. In all cases 
the deformed sample becomes doubly refracting; this property disappears after 
heating without pressure. When the deformed sample is interposed between the 
polarizers and illuminated by a parallel light beam there appears in the field of vision 
a dark cross, whose position does not change when the sample is rotated about the 
axis perpendicular to its plane (Fig. 1). This is a typical picture of a spherolytical 


Fig. 1. A polyvinylchloride plate between crossed polarizers. On the Fig. 1b the sample is rotated 45 0° 


structure; it means, that the principal optical axes of the various parts of the sample are | 
oriented radially and tangentially, and that in one of these directions there exists some _ 
orientation of the molecules. This effect is obviously connected with the „flow“ ae 
the chain segments described above. 

By means of a quarter-wave plate it was possible to estimate the difference An 
of the refractive indices in both principal directions by measuring the phase difference 
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Ag of the elliptically polarized light. It is known, that 


where À — the wavelength of light used (assumed to be 5900 A). 
d — the thickness of the sample. 
The calculated value An varies from zero at the center of the sample to 10-4 
near the edge. 


The average refractive index of a sample was also estimated by the measurement 
of the Brewster angle; it was found to be 1.53. 


c. Electrical charges 


The samples of polymer after forming under the pressure obtained a surfece charge. 
This charge was measured by means of the lifted electrode method described by 
Eguchi (1925). This method gives an indication of the average surface charge density 
on both sides of the sample, and the total electrical charge. 

It was found in this way that in all cases an assymetry of surface charge exists, 
which is influenced by the body which is in contact with the sample during the formation 
process. In the case of the aluminium or tin foil (case A) the total charge of the sample 
is negative. At first both sides of the sample have negative charges of the order 
200 CGSE, which diminishes with time but not in the same manner. The negative 
charge on one side decreases slowly, with a time constant of a few months. The nega- 
tive charges on the other side decrease, at first rapidly to a small value, or change 
sign to positive, and then decay slowly. ہ‎ 

In the case of the bakelite or copper plates (case B) the total charge is nearly 
zero; the charges on both sides have the opposite sign and are of comparable value. 
These charges also decrease with a time constant of the order of a few months. 

A more detailed picture is obtained by the measurement of the surface charge 
distribution. The first observations made with the aid of an edge counter used as an_ 
electrometer as described by Sujak (1957), showed that the charge is not uniformly 
distributed on the surface. This method as yet cannot give any quantitative results; 
the “lifted electrode“ method was more convenient. The lifted electrode was made 
smaller (0,5 cm?) with a guard ring lying ont he sample (Fig. 2). By changing the 
position of the electrode and the ring on the sample one can measure the charge 
distribution. The results are represented in form of “map“; a typical one is illustrated 
in Fig. 3. l 

Analysis of these „maps“ permits one to make several conclusions: 

a) On one side there exist areas having charges of different sign. 

b) The maximum of the charge density lies at the boundary of the sample. The 
negative charges may, also have their maximum at the center. 
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c) In case A one observes at first uniformly distributed negative charges on both 
sides. A few hours later there appear positively charged spaces. 

d) The spaces charged positively at one side are accompanied by a negative 
charge on the opposite side. 

This latter conclusion is an important one-it may be assumed that in the polymer 
there exist spaces having a semi-permanent electrical moment, as well as electrets. 
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Fig. 2. The device for the measurement of the charge density distribution. A-upper electrode, B-guard 
ring, C-lifted electrode, D-electrometer. 


The picture of the charge distribution does not change with the time (excluding 
case A in the first period), but the charges decrease. The maximum of the charge 
density observed was 27 CGSE, which is greater then the value observed in electrets 
(Gubkin, Skanavi 1957). The samples were preserved in desiccator, wrapped in foil. 
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Fig. 3. The distribution of the charge density on one polymer sample. The density given in CGSE units. 
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§ 3. Discussion of results 


The observed effects: deformation of the sample, double refraction and electrical 
moment can be correlated with one another. 

According Skinner, surface charges on a polymer in contact with another body 
increases owing to the passage of the charge carriers through the surface of the contact. 
One would expect that in our case (identical material on both sides of the sample) 
both surfaces would have charges of the same sign and of the comparable value such 
as at the beginning of the measurements in case A. The assymetry of charges can be 
explained if we assume the existence of a dipole orientation in the clamped polymer. 

On the basis of X-ray investigations Mark (1943) concludes that in ordered regions 
the dipole groups of the polymer and the chain axes have a parallel orientation with 
the dipole axes perpendicular to the chain axis. Such an orientation seems to appear 
in the conditions of our experiment. The thermal agitation of segments and the 
existing pressure tends to align the chain segments in a more ordered position. 

Following Gurnee (1954) we assume such parallel oriented chains as having 
100% orientation. If we assume the dipol moment of a monomer as roughly u = 2 
Debye units, then one can calculate the degree of the dipole orientation 


LEE OO à 
O dip = Pee == Nu = 107 
where P is the existing electrical moment of volume unit, equal measured surface 
charge density o. i 
= Ny is the calculated moment in the ordered region N is the number 
of the dipoles in the unit volume. 
In a similar way we can evaluate the degree of the optical orientation by calculating 
the difference An = mı — ny of the refractive indices parallel and perpendicular to 
the chain axis. They are related to polarizabilities of the chain segments by the equation 


(Gurnee 1954): 


E 


max 


27 
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where n = mtr. is the average refractive index. 


Polarizabilities a and @’, can be calculated by taking the data of تبیہ ادن‎ 
of bonds (Volkenshtein 1955). Assuming, that the tetrahedral angle is roughly 120 
one obtains for the’ static model of the completely oriented polymer: 


ay — aj; =3 x 107%, Anma = 0.35 

If we assume, that the H—C—H and H—C—Cl groups rotate about the axis 
perpendicular to the chain axis, then one obtains 
f ay — e = 2.7 X 10-24, Anma = 0.31 
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Dividing the observed value An by the calculated one, we obtain the degree of 


the optical orientation 


An 


Te‏ == وہ0 
max‏ 


~ 0.3 X 10-3, 


which is comparable to the degree of dipole orientation. 

The radial distribution of the charge, as mentioned above, is such, that they 
are larger at the edge of the sample; the observed double refraction varies in a similar 
way. This effect may be connected with larger flow velocity at the edge of the sample. 

In the case A we must assume that there exist charges of both kinds: 1) the contact 
charges (negative) which diminish rapidly, and 2) the polarization charges, which 
decay slowly. 

The proposed mechanism for production of the static charges should be confirmed 
by X-ray investigation and dielectric constant measurements, and also by examination 
of other polar and non-polar polymers. 

I ame greatly indebted to Professor B. Makiej for the support he has given to my 
investigation. The author also wishes to thank Professor H. Kuczynski for the polymer 
sample, and Professor B. Sujak for help in measurement on his apparatus and for 
many helpful discussions. 
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The phosphorescence decay curve of acridine yellow in gelatin after long preexcitation 
in linearly polarized light has been investigated. It was found that the decay constant of 
the luminescent group of centres increases exponentially with the preexcitation time. 
An increase in the decay constants conditions the decreases in the total light of phosphores- 
cence in the deformed (preexcited) state. These changes are caused by the appearance, 
during strong preexcitation, of unstable isomers of the dye. 


1. Introduction 


_ The object of this paper is to investigate the influence of preexcitation on the 
phosphorescence decay curve of acridine yellow in gelatin. 

As a result of long illumination in intense light a number of organophosphors 
change their properties. Such a state of the phosphor will be called the deformed state 
and will be denoted by P. The state of the same phosphor not subjected to preexci- 
tation will be called the undeformed state and denoted by ۰ 

Frölich and Szalay (1948) distinguished two characteristics of the deformed state 
according to the ratio of the intensity of the emitted phosphorescence in the deformed 
state to the intensity of the phosphorescence in the undeformed state for the same 
conditions of excitation. If the intensity of phosphorescence of the state P is greater 
than the undeformed state P®, then we speak of a positive P+ deformed state and 
in the opposite case, of an negative P- deformed state. 

It turns out that the character of the deformed state of the phosphor can be changed 
from P+ to P- by a change in the temperature. A measure of the deformed state 
(according to the phosphorescence intensity) is the ratio I(P)/I(P°) = A where I(P) 
denotes the intensity of phosphorescence in the deformed state P and [(P°) is the 
intensity in the undeformed state. The experiments of Frölich and Gyalai (1936) 
showed that the state P and A change with the time ellapsed from the moment the 
preexcitation is interrupted. The deformed state P tends to the undeformed state 


f P® and A> I. 
(93) 
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The numerical value of 4 depends on the time and intensity of preexcitation, and 
on the direction of polarization of the preexciting and exciting light. For some dyes 
the value of A is greater if the electric vectors of the preexciting and exciting light 
are parallel than if they are perpendicular to each other. For other dyes, the opposite 
is the case. The value of A depends also on the dye concentration. 

The investigations of Lewis and his coworkers (1941, 1942, 1943) showed that 
there exist molecules: 

a) whose absorption optical axes in the deformed state (preexcitation is obtained 
by excitation of the phosphor in linearly polarized light at the maximum of the 
absorption band are parallel to the electric vector of the preexciting light, and 
such molecules 

b) whose absorption optical axes in the state P (with the same preexcitation as 
before) are parallel to the electric vector of the linearly polarized preexciting light. 

In case (a) we obtain maximum absorption of linearly polarized light if the electric 
vectors of the preexciting and absorbed light are parallel. In case (b) the absorption 
maximum is obtained if these vectors are perpendicular to each other. The above 
described properties induced Lewis to divide the molecules into two groups. In the 
first group he included molecules whose properties are described under point (a) and 
he considered the orientation of their axes to be normal; in the second group he 
included molecules with properties described under point (b) and he considered the 
orientation of their axes to be anomolous. 

On the basis of the division introduced by Lewis, we can distinguish the deformed 
state P with normal or anomolous orientations of the absorption axis of the molecules. 

The deformed state is maintained for a longer time after interruption of the 
preexcitation. In gelatin, for example, for Rhodulin-Orange-N the deformed state 
disappears after 1 to 3 hours depending on the dye concentration. In rigid solutian, 
for some dyes the orientation of the molecules is stable and can be destroyed only 
by dissolving. 

Fröhlich and Szalay (1948) attempted to give a full explanation for the deformed 
states P+ and P- by assuming that during preexcitation the molecules become oriented 
parallelly or perpendicularly to the direction of the electric vector of the linearly pola- 
rized preexcitation light. This assumption leads to some conclusions which the authors 
have not yet been able to establish experimentally. This is the reason for undertaking 
investigations of the deformed state. 


2. Experimental Procedure 


The investigation of the deformed state was made on acridine yellow in gelatin. 
The phosphor was prepared in the following way: An unexposed film strip (Ultra- 
Rapid) 1.80 m in length was fixed and dried. It was then placed in a bath of an aqueous 
solution of acridine yellow, rinsed in running water and dried in a dark room. The 
concentration of the dye in gelatin was determined by means of a Pulfrich colorimeter. 
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A linear phosphoroscope was used to investigate the phosphorescence deca 
curves. A detailed description of the method of measurement of the decay curves “ 
of the linear phosphoroscope is given the paper of Frackowiak (1957a). In order to 
increase the intensity of preexcitation an additional light source and optical arrangement 
was connected to the phosphoroscope. A diagram of the set-up is shown in Fig. 1. 


See 
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Fig. 1. Diagram of phosphoroscope. Z — light source, L — lenses, W — water filters, P — polaroids, 
K — trap, V —voltmeter, D — diaphragm, C — photographic plate, M — motor, T — tachometer, 
F — phosphor. 


; 
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The phosphor is preexcited by both light sources A and B (Fig. 1). The measure- 
ments of the decay curves are made by means of arrangement A. The phosphoroscope 
speed during preexcitation is considerably smaller than during the measurements of the 
decay curves. In such an arrangement oneobtains interrupted preexcitation. Since, 
however, the duration of the deformed state is of the order of hours, and the modulation 
frequency of the light is of the order of sec”, it can be assumed that the deformed 
state during the time between two successive illuminations of the phosphor does not 
effectively change. This permits the use of the preexcitation time + just as in the case 


of continuous preexcitation. 
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t can be calculated from the phosphor speed and geometry of the arrangement. 

Henceforth in this paper we shall use only the time t calculated in this way. 

After a given time of preexcitation, the preexcitation is interrupted and the speed 
of the phosphor is increased. Thirty seconds after the interruption of the preexcitation 
source A is switched on and the decay is measured. The measurements were made 
under the following conditions: 

1. Room temperature 20°C. 

2. The phosphoroscope speed during preexcitation was (382 + 60) cm/sec. 

3. The phosphoroscope speed during excitation was (1150+ 60) cm/sec. 

4. The time of exposure was 60 sec. 

5. The time between the interruption of the excitation and the moment of obser- 
vation of the decay curve was 5 X 1074 sec. 

6. The preexcitation times were t = 0, 1, 3, and 6 min. 

7. The preexcitation and excitation light was white, linearly polarized (with 
parallelly oriented electric vectors). 

8. The intensity of the preexcitation light was twice as large as the excitation 
intensity. 

9. The measurements of the decay curves for different times of preexcitation 
were made for one and the same phosphor. 

10. Three hours passed between two successive preexcitations, the time of the 
decay curve measurements being included. 

11. The time of the decay curve measurements for each preexcitation time t 
did not exceed 5 min. 

12. The excitation time was tọ = 1074 sec. 


3. Experimental results 


The components of the decay curve of phosphorescence parallel 7|, (t) and perpen- 
dicular 7 , () to the electric vector of the linearly polarized excitation light was expanded 
by a well-known method (e. g. Frackowiak 1957a) into a series of exponential functions.. 
The components of the decay curves and their exponential functions are given in Figs. 2, 
8, 4, and 5. 

The decay constants determined from these curves for the individual liminescent 
groups of centres and their limiting degree of polarization are collected in Table I. 

Knowing the components 7, (£) and I, (f) one may determine the decay curve 
of phosphorescence I(t) emitted in all directions from the relation 


I) = 1,0 + 21,0 8 
The decay curve I(t) found in this way can be compared with the theoretical 


curve for isotropic rigid solutions (Jabłoński 1954, 1955, 1957) for the simplified 
model of the centre 


2 k—1 
1) = BY ساس‎ ex) 7a (2) 
: 1 | 
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Table I 


Acridine yellow in gelatin. Dependence of the decay constants of groups of centres and their limiting 
degrees of polarization on the preexcitation time 


Preexcita- Decay In decay Limiting degree 
tion time constants constants of polarization 
T, sec 1 nl, Pk % 
0 5.86 X 103 8.67 36.9 
6.50 x 10? 6.48 25.7 
1.60 x 10? 5.07 12.4 
سی رت‎ |i 
7.60 x 103 8.93 21.2 
2 1.05 x 103 6.96 31.7 
1.80 x 10? 5.19 10.5 
6.95 x 108 - 5 269 
2 1.10 x 108 7.00 12.7 
3 1.90 x 10? 5.25 8.5 
6 1 | 77x10 | 8.96 31.8 
2 1.42 x 103 7.26 28.3 
3 2.50 x 10? 5.52 8.9 
In J 


1 2 3 4 5 6x10 ser t 


Fig. 2. Acridine yellow in gelatin. Logarithm of phosphorescence intensity illumination time tT = 0 min 
—o— parallel componenet and its exponential functions, 


—X— perpendicular component and its exponential functions. 


1 Measurement made under somewhat varied conditions. 
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1 2 3 4 5 6x10 sec t 


Fig. 3 type same as above, but 1min. 


where B — is a proportionality coefficient 
Lait ا‎ EF 
y, — is the probability of an electron transition accompanied with emission 
of photoluminescence in the k-th group of centres 


W, — is the probability of internal decay and decay by molecules of the solvent 


W, — is the probability of decay the unexcited luminescent molecules in a centre 
of the k-th group 


y=n-v 
n —is the number of luminescent molecules in 1 cm? of solution 


v— is the volume of the active sphere of the centre 
t— is the time 


k —is the index of the kind of group to which the given centre belongs. 


1 2 3 4 5 6x10 sec t 
Fig. 4 type same as above, but 3 min. 
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n F. 


1 2 3 4 5 6x10 sec t 
Fig. 5 type same as above, but 6 min. 


If it is assumed that there is no decay in the solution under investigation, then, 
I, = y, Eq. (2) passes over into (3) 
yk-1 
I(t) = B 2 yg (1 — e kt) — e” (3) 
= (k — 1)! 

All the constants in the theoretical formular (3) except for v can be determined 
from the analysis of the experimental decay curves. The value of » can be found by 
a trial method. | 

A value of » is chosen so that the theoretical curve best coincides with the experi- 
mental curve. 

Fig. 6 shows the experimental and theoretical (3) decay curves of phosphorescence 
of acridine yellow in gelatin for different preexcitation times T. 

The values of v found in the above-described manner and the radius of the active 
sphere calculated from the relation v = nv are given in Table II. It is worth noting that 
the mean value of the radius of the active sphere is in agreement with the value obtained 
in the previous investigations of the decay curves (Frackowiak 1957b). 


4. Discussion 


All decay curves, as may be seen from Figs. 2, 3, 4, and 5, can be decomposed 
into three exponential functions irregardless of the preexcitation time. The theoretical 
curves shown in Fig. 6 have (within the limits of error) the same values v and the 
same values of k (k = 1, 2, 3). These curves, for different preexcitation times, differ 
only in the proportionality coefficients and decay constants. It thus follows that: 

1. both in the deformed and undeformed states the same groups of centres 
(groups of centres of the same population and the same active sphere radius) take part 


in the process of phosphorescence; 
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Table II 


Acridine yellow in gelatin. Variation of », radius of active sphere, and totel light of phosphorescence 
with preexcitation time 


Preexcita- 2 Radius of active : 
à : Light sum 
tion time sphere 
— | J 2 ٦ T 
T min v r(À) ت7‎ A 
k 
0 7.9 221 3.02 
l 7.9 221 2.47 
3 | 10.0 23.5 2.84 
6 | 8.0 22.3 کا ا‎ 
ee 0 28 + 
Mean value 8.45 22.5 
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200 


700 
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2. the decay of phosphorescence in the deformed and undeformed states is des- 
cribed by Eq. (3), differing only in the proportionality coefficient and decay 
constants. 

As may be seen from Fig. 6, the intensity of phosphorescence decreases with an 
increase in the preexcitation time. An exception is the curve denoted by the number 3. 
The measurement of this curve was made under somewhat different conditions 
(24 hours after the measurements of the other curves). We thus ascribe less weight 


0 7 3 2 T min 


Fig. 7. Acridine yellow in gelatin. Logarithm of decay constants of groups of centres vs preexcitation 
time 2 


to this curve. We may therefore state that, in the case of our experiments, acridine 
yellow in gelatin changes into the deformed state P- as a result of preexcitation (in 


accordance with the definition of Frölich). rs 
In analyzing Table I it may be noted that with a change in the preexcitation time 


limiting degree of polarization P, of a given group of centres and their decay constants 


I, also change. ۱ K ہیوت‎ 
: The limiting degree of polarization changes irregularly so that it is difficult to draw 
conclusions from these changes. On the other hand the decay constants T, change 


very regularly with the preexcitation time. : 
Fig. 7 show the changes of the natural logarithms of the decay constants I, as 


a function of the time of preexcitation 7. 
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As seen from this figure, these changes are represented by a series of parallel 
straight lines. These lines can be described (in the limits of the preexcitation time 
investigated) by the formula 


I, = y, (4) 


where I, — is the probability of electron transition accompanied by the emission 
of phosphorescence in the deformed state, 
y, — is the probability of electron transition accompanied by the emission 
of phosphorescence in the undeformed state, 
æ — is a constant which we shall interpret later, 
t — is the time of preexcitation. 
The observed changes of the decay constants of a group of centres can be employed 
(analogously to the characteristic of the deformed state given by Fröhlich (1948) to 
characterize the deformed state. 

If Ty, > 1 or I,/y, < 1 we shall speak of a positive or negative deformed state 
with respect to the decay constants of a group of centres. We denote it therefore 
by (+ or —). 

The observed process of the growth of the deformed state, as seen from Eq. (4), 
is exponential. Fröhlich and Gyulai (1936) observed the extinction of the deformed 
state P* with time after the interruption of the preexcitation. He also noted the exponen- 
tial character of this decay. 

The decay curves of phosphorescence of the undeformed state P® and the deformed 
state P+,- are described by Eq. (3) and differ, as already mentioned, only in the 
proportionality coefficient and in the decay constants. The increase of the decay 
constants J’, in the state ۶٣۴ with a simultaneous decrease in the intensity of the 
emitted phosphorescence indicates the existence (in this state) of a “quenching“. 

Because of this, let us ecamine the assumptions on which basis Jabłoński (1957) 
derived Eqs. (2) and (3). The most important assumptions from the point of view of 
our discussion are: 

1. Eq. (3) is valid only in case there is no phosphorescence decay, that is, if 
W;= W, = 0. In this case the total light of phosphorescence should be constant 
if the absorption or the ratio of the fluorescence to phosphorescence efficiencies does 
not change. 

2. The probability of transition from a fluorescent level F to a metastable level M 
is proportional to the transition from the level M > F and equal to By,. 

As may be seen from Fig. 6 (curve 1), the theoretical curve described by the 
formula fits very well to the experimental curve of phosphorescence decay in the 
underformed state P®. 

Thus if there exists a “quenching“ then it appears only in the undeformed state 
P--*, It thus follows that the decay curve in the deformed state must be described 
by Eq. (2) where پ7‎ = y, + Pp, and where gk denotes the سو و‎ proportional 
to time of any process quenching the phosphorescence. 
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Making use of Eq. (4) and (2), we find the formula for the phosphorescence 
decay curve in the undeformed state P+ as a function of the preexcitation time 


Ts 
7 = ہے 7چ2-‎ pay ار‎ 
| (t) = Be 3 7701 یی ا‎ td (5) 


As may be seen from Eq. (5), for a fixed preexcitation time + the phosphorescence 
decay curve in the deformed state ۶-۳۰ is described by a formula of the same as (3) 
but differing from it only in the proportionality coefficient and the decay constants. 
Other decay curves shown in Fig. 6 (curves denoted by the numbers 2, 3, and A) 
are compared with the theoretical curve (5). 

In order to interpret the constant œ occuring in the formula we shall investigate 
the total light of phosphorescence L(t) in the state ۰م‎ as a function of the prexcitation 
time t From (5) we obtain 

O0 


gor | OR ٤ 


ہے 
k k‏ 


0 


where J, (0) is the intensity of the phosphorescence emitted by the k-th group of centres 


at the time the excitation was interrupted. From the experimentally obtained decay 
Le (0) 


curves we may determine J, (0) and T, or yẹ. The numerical values >} 7, are 
k k ! 


given in the last column of Table II. ۱ 
In Fig. 8 are given the natural logarithms of the total light of phosphorescence as 
a function of the preexcitation time. As seen from Fig. 8, the function L(t) decreases 
exponentially with an increase in 7. From the slope of the straight line shown in Fig. 8 
we can determine the decay constant Û of the function L(t). The absolute value of this 
constant is B = 16.0 x 10-4 sec. 
It follows from Eq. (6) that the function L(t) is not a simple exponential function 
because of the dependence of J’, on t defined by Eq. (4). In our experiments, however, 


inl 


6 Tmin 


Fig. 8. Acridine yellow in’ gelatin. Logarithm of light sum of phosphorescence vs preexcitation time, 
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the value of the expression e~7#* appearing in Eq. (6) is contained, as has been calculated, 

in the numerical interval 0.4 < ۸م‎ < 0.5 for all times t and decay constants I. 

This factor thus has such slight effect on the function L(t) that in the interval of 

preexcitation time Wwe are investigating, it can be treated as an exponential function. 

In a similar way the constant « occurring in Eq. (4) can be determined from the slope 

of the lines shown in Fig. 7. The value of this constant is œ = 8.3 x 10-4 sec. 
Comparison of the values of both constants ( and f) gives 


B = 2a (7) 


which proves the correctness of Eq. (5), and (6). 

It thus follows that we can interpret the constant œ as a quantity proportional to 
the transition from the undeformed state P? to the deformed state ۶۶۰۶۰۰ Let us consider 
at present the energetic scheme of the measurements of phosphorescent molecules 
given by Jablonski (1935). Let N, M, F denote the normal, metastable, and fluorescent 
levels in state P® and N’, ۸7, F’ the corresponding levels in state P~*. 

In the undeformed state, at room temperature, the probability of an electron 
transition accompanied by the emission of phosphorescence is practically equal to 

` the transition M —> F and equal in our case to y. (It is sufficient to discuss one lumines- 
cent centre of any group and we omit the index k). Similarly, in the deformed state 
the probability of transition accompanied by the emission of phosphorescence is 
equal to the transition M'— F’ and, by Eq. (5), equal to I. 

In the undeformed state P? the probability of the F> M transition is, as already 
meantioned, equal to By. On the other hand, in state P+, by (5) and (4) the probability 
of the transition F’ -< M” is equal to Be ®” [ = Be“ y. Comparing the last relation 
with Eq. (4) we see that in state P+ the probability of the transition M’ > F’ decreases 
with the increase in the probability of the transitions F’—> M” (according to the 
inverse of the law (4)). Owing to these changes in the state P™* the phosphorescence 
efficiency decreases and the fluorescence efficiency increases simultaneously and, 
as a result, there is observed a decrease in the total light of phosphorescence of state 
7۲۳۰ with respect to P®. 

The phenomenon appears to proceed as if in state Pt there exists a „quenching“, 
which, as shown above, really does not exist. 

It may be noted that this phenomena does not depend on whether the absorption 
in the deformed state changes in relation to the undeformed state, or whether it remains 
the same. It thus follows that the changes in absorption which can take place in the 
deformed state are a characteristic feature independent of case in the intensity of 
the emitted phosphorescence. S 

Comparing Eq. (4) and (6), we note that in the case of the decay constants of 
a group of centres in the deformed state decreasing according to the inverse law in 
relation to (4) we obtain an increase in the total light of phosphorescence. In accordance 
with the nomenclature adopted by us we would obtain in this case the deformed 
state P+, It thus follows that the second sign describing the character of the deformed 
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state, owing to changes of the probability of the transition M’—> F ', defines the first 
sign. It may therefore be assumed that the deciding factor conditioning changes of 
the inensity of phosphorescence in the deformed state are ‘changes in the probabilities 
of the transitions M’—> F’ and F' > M. 

In order to explain the occurrence of the deformed state it may be assumed, as 
one of many possibilities, that during strong and long preexcitation there take place 
certain reversible structural changés (unstable isomers) in the luminescent. molecules. 
The centres formed from undeformed molecules as well as from deformed molecules 
(isomers) therefore take part in the process of phosphorescence in the deformed 
state. | 

If the deformed molecules (isomers) differ in certain characteristics from the 
undeformed molecules, then the deformed state should also differ from the 
undeformed state. i igh? 

Since, as already mentioned, after preexcitation both types of centres: (deformed 
and undeformed) take part in the phosphorescence, then the properties of the deformed 
state are the result of the superpostion of effects, and may more or less differ from the 
properties of the undefomed state according to the number of deformed centres. 

As a result, for a given intensity of preexcitation and preexcitation time it should 
be expected that the effect of the deformed state depends on the concentration of the 
luminescent molecules. The investigations of Frélich and Gyulai (1957) showed ithe 
existence of such a dependence. 

It seems necessary to discuss in detail the system composed of the undeformed 
and deformed molecules and to make further investigations of the various properties 
of the deformed state. 


6. Summary 


The décay curves of acridine yellow in gelatin have been investigated as a function 
of the preexcitation time. From the analysis of the decay curves the decay constants 
of the individual groups of centres and their limiting degree of polarization has-to be . 
determined. Changes in the decay constants of groups of centres the preexcitation 
time has been expressed in a formula as a function of time. By means of this formula 
and the formula, of Jabłoński for the decay of phosphorescence in the undeformed 
state a formula for the change of the total light phosphorescence in the deformed state 
has been given as a function of the preexcitation time. It has been shown by comparison 
there is agreement between the theoretical and experimental decay curves of phospho- 
rescence in the deformed and undeformed states. The radius of the active sphere 
of centres in both states has been calculated. The radii of the active spheres in both 
states prove to be the same. A new characteristic of the deformed state dependent 
on the changes of the decay constants and total light of phosphorescence has been 
introduced. It is shown that the reduction of the total light of phosphorescence in the 
deformed state is connected with the increase in the ratio of fluorescence to phospho- 

j i 


\ 
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rescence efficiencies. This increase arises as a result of the increase in the probability 
of transition from the metastable to the fluorescent state with a simultaneous reduction 
in the probability of transition from the fluorescent to the metastable state. 

The occurrence of the deformed state is explained by the appearance in the 
phosphor, as a result of preexcitation, of unstable luminescent isomers of the dye 


under investigation. 
The authors sincerely thank Professor A. Jablonski for his valuable discussions 


and comments in the course of this work. 
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The scheme presented here permits calculations of the influence which the positive 
centres distributed along a potential box of the linear free-electron model have on the 
energy states of this model. The numerical calculations included the configurational 
interaction for the two and four centres problems (ethylene and butadiene in the z-electron 
approximation). The computations were done by assuming an equidistant distribution of 
the centres in the model and that the box radius tends to zero. The results obtained for the 
case of the lowest excited singulet states are in surprising agreement with the correspon- 
ding states calculated by means of the LCAO MO model. At the same time it was found 
necessary to assume a very strong screening of the nuclei. This may be considered as, 
a consequence of the geometry of the problem. : 

The scheme was applied also to calculation óf the positions of the absorption maxima 
corresponding to the lowest electron transitions in Simpson symmetrical cyanine dyes. 
The results remain in especialiy good agreement with experiment if it is assumed that the 
nitrogen atom localizes the free-electron pair. 

It is pointed out that if one includes the effect of the interaction of nuclear charges 
on the electrons of the model, it is found that the energy states group themselves in zones 
analogous to the energy zones known from the band theory of solids. 


The purpose of this paper is to calculate the influence of the nuclear charges 
apon the energy states of the free-electron model. We assume that the electrons are 
moving in a potential box extended along the free-electron bond path, interacting with 
one another and remaining in the field of the positive point charges Ze. These charges 
represent nuclear charges of the atoms of a conjugated chain obtained by removing 
the electrons from the atomic cores; these electrons are further considered as free- 
electrons. Thus, in the case of the carbon atoms, for instance, if we assume that each 
atom contributes only one -electron to the electron gas of the molecule, we have 
Z = 1; we take Z = 4 if it is assumed that the ø as well as the x electrons form the 
electron gas in the potential box of the model. (For details of the or — o treatment in 
the free-electron model see the previous work of the author (1957b)). 

In this paper we present a scheme which permits calculations of the energy 
states as a function of the magnitude of Z and of the distribution of the positive centres 
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of the model. The numerical calculations for several simple particular cases are also 
given. 

In our further considerations the term “free-electron model“ will be used, although 
it corresponds. only to the treatment of the problem in the zero-order approximation. 
We shall presume that the picture in our present scheme differs from this original 
treatment in that it takes into account the interaction effects which we have mentioned 
above. 


§ 1. Computational scheme 


The Hamiltonian operator of our problem now has the form: 
H =H? +H! + A (1) 


; CBA h2 : 
where H° is the kinetic energy operator equal to — >) A 72; the next term H’ is the 


Coulomb energy operator of the interaction’ between the electrons and the nuclei: 


ieee za 3 =r (2) 


where r, indicates the position of uth electron, &, — of jth nucleus, and finally H” 


is the mutual electron interaction operator equal to >) >) پت‎ (u < (۰ 
v K 

We restrict ourselves to the following Apa as regards the geometry of 
the problem: the electrons move inside a cylindrical potential box; the atomic nuclei 
(usually the carbon nuclei) are situated along the box axis at equal distances from 
each other. The box corresponds to a straightened chain of the linear molecule with 
a conjugated double-bond system. The length of the box exceeds the length of the 
chain by one interatomic, distance at both eneds (Fig. 1). 


Fig. 1 The butadiene box 


Each of these assumptions has already been discussed in detail elsewhere. It seems 
unnecessary to repeat this discussion now. i 
The molecular orbitals of our problem are the eigenfunctions of the operator H°. 
As has been shown previously their form, if the boundary conditions are taken into 
account, is: i 
na 


P : 
Qn, (Y) = 8 sin Fre" N; J (wr) cos 18 (3) 


r =r (z,r, 0) 
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where ! = 0 corresponds to the o, and 1 = 1 — to the z-electrons. The energy 
values furnished by H° are the the well-known eigenenergies of the simple free- 
-electron model. 

We further regard HT and H” as the perturbation operators. We compute the 
energies corresponding to these components of the Hamiltonian operator by means 
of @,, taking as the function in the zero-order approximation for each of the states 
suitably antisymmetrized products of g. 

The effect of H™ on the energy states of the free-electron model has been discussed 
previously (see Olszewski 1955, 1957b — henceforth denoted as I, III). Let us now 
consider the contributions introduced by H7. In general, these are sums of integrals 
of the type 


core Ze? 

CSS E n, (r ae er کا‎ LR dt 
D frere a ۵ 

= — Ze? > | pr) 4nG (r, a) qur (F) dr (4) 
7 


G is the Green function of the three-dimensional Laplace-Poisson equation. If we 
assume that j th nucleus which we are now considering lies on the box axis at a distance 
a; from the origin of the free-electron path, i. e. aj, = |a| = dj, dj, = 0, then the 
integral representation of G in cylindrical coordinates takes the form: 


+00 ہہ‎ 2r 
1 ix|z—a; | 1 : 
G(r, a) = = / rte finena | oe- oa (5) 
—00 0 0 


The last integral in (5) is equal to unity and may now be omitted. Inserting G into (4) 
we may calculate the integral which is of interst to us. From the integration over @ in (4) 
we see immediately that [°° vanishes unless J = J’. 

The elementary trigonometrical transformations of the factors depending on z 
in the molecular orbitals product under the integral sign (4) allow us to represent each 
of the j components of 1% in the form of the difference of two integrals of the type 


L R ہو‎ : 
D g= 1 f il ie cos cz ° N? Ji (wr) cos? ۰ 4nG(r, aj) dz rdr dð )6 
02 | 


We calculate (6) for the general case c 0, and for c = 0. 


We integrate first over z, next over x, k. In order to show that such a treatment is correct we shall 
prove that the integral | 
Me 


اچ سن 
abl ane Mee eee (a)‏ 
gh = KS Ab‏ 
co 0‏ — 
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is uniformly convergent in the whole region of z (O, L). Indeed, carrying out integration over x we obtain: 
O0 
i |زہ-ھ | طس‎ Jy (kr) dk (2’) 
0 


co 


Without further computations it may be seen that, in view of the convergence of f Jo (kr) dk and owing 
0 


to the fact that for any arbitrary z the function f (k, z) = e~*! ز“-ھ‎ | diminishes monotonically, and for 
مز‎ = 0 is smaller than a constant C independent of z, the conditions of Abel’s criterion for the uniform 
convergence of the integral (a, a’) are fulfilled. 


We obtain the following results for S: (case c ¢ 0) 


R 
Ze? ; 
S = Te { [sin ca; {Si (cay) F Si [c (L — a;)]} + cos ca; {Ci (cag) + 
. i 
+ Ci [c (L —a;)} | Ip (cr) + 2 cos ca Ko (cr) — Pa; (cr) — ہو-رط‎ (er). (7) 


N J? (u,r) rdr, 
forsc =) 0: 


R 
2 : 
S= — < fare sh (2 i) gee sh pen Ni J? (ur) r dr. (8) 
0 


In the computations we took into account the boundary conditions which the free- 
-electron molecular orbitals fulfill. The upper sign in (7) holds when cL is an even 
multiple of x, the lower — when it is an uneven multiplet. 

P, (cr) has the following meaning: 


1 ۱۶ء / de mil!‏ ری 
)9( ےہ ت0 Gi) = Nae D (— 1) (cay‏ 


and K, (cr) has the convenient expansion: 


2m‏ ا 
K o (cr) = — log = æ) I, (cr) + > — y(m + 1), (10)‏ 


1 Nj are the normalization factors of the J, (u, r) functions. Their values are 
2.724 3.514 
0 ae ; 
R Ê R & 
R is the box radius. 
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where 


1 


y, y — Eule’rs cosnstant.‏ — ا 


1 
کا ہے ا اچ با سا 


For large c, P (cr) tends quickly to — SL [Z (er) -- 1]. 
a? c 


Let us now consider the case where the ratio of the radius of the box to its length 
(or, what is equivalent for constant L, the radius alone) tends to zero. It is generally 
accepted that this assumption in the calculations for a given molecule is the more 
sound, the larger the length of the molecule. It seems to be of particular importance 
to emphasize that the results thus obtained on the basis of the above assumption are 
especially valid in the case of very long free-electron paths. 

To simplify our considerations we may assume that the charge inside the box 
is uniformly distributed around the axis, i.e. instead of J? (u, r) in (7), (8) we have 
unity. By means of this assumption calculations of the mutual interaction energy in 
several conjugated systems have been performed (see I, III). 

The 7%" integral representing the energy of interaction of the electron in the 
nth state with the nuclei of the molecule then takes form: 

1 


J. 


E {ci |a — à) آ2‎ Jaa (252 « cos =e a 


E 2771 1 : 70 | Inn 
= 2.008 2 ay flog? +} — lim og R — 7} + 200 7 a; log ہے‎ | 


ENS Ez 4 log (L — a) — 2 lim log R + 1+ 2log2 + 
R-0 


: 2an r 
The above expression tends to infinity, unless D ( 1 — cos 7 à) = 0, The ful- 


filment of this condition is, however, not indispensable for the use of equation (11). 
In most cases we are interested in the difference of energies of two states. Therefore, 
in order to make possible the computation of the difference of electrons-nuclei 
interaction energies in the state n and m, respectively, when R/Z > 0, the relation 


2 27 71 
D E = Ope COS ate «| =Ù. (12) 


7 
the contribution to the computed excitation energies, is then 


should occur. Thus, تس‎ : 
log m. We also assume, as is essential 


given only by the expressions with Si, Ci, log n, 
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in the case of a very thin tube, that the configuration of the nuclei does not change 
during the excitation. 

We compute the integrals of the type If (n m), in a way similar to that 
of IS9", Such integrals appear in cases in which the configurational interaction is 
included in the calculations and take the form 


1 
lim ۸3 z (Seen A Speminem | = 


R+0 S 
= بہت‎ sn (m — 7) La fs [em -« Ta] Fsi Cr -a + 


9 


+ cos (m — n) Z ylei ۲ — n) = 7 رت‎ ۲ — n) > (L -a |+ 
sin (m +n) = رہ‎ ‘si | +n) + 7 AS | +n) = L- a) + (13) 
+ cos (m + nZ alci (m +n) ا‎ i ۲ +n) > &a) |} 7 
+ 2 cos | — n) + | log F (m— n) +2 7 | +n) Ta | log + (m + n) | ; 


The upper signs correspond to even m + n, the lower — to uneven m + n. The integral 
(13) only has sense if a condition analogous to (12) is fulfilled: 


mr SDT 

2 sin TT G'sin—- aj = 0 (14) 

j 
The odd m + n appear only in the case of a linear free-electron model where the 
configurations of different symmetry (4 and B) interact. As is known from group- 
-theoretical cosiderations, the integrals caused by this interaction should vanish. 

The integrals representing the influence of the nuclei are inversely proportional 

to L, as were the electron interaction integrals. 


LE $ 2. Simple applications 


Let us apply the scheme presented above to some simple cases such as ethylene, 
butadiene and several symmetrical cyanine dyes investigated by Simpson (1948; 
McGlynn and Simpson 1958). We shall assume in the calculations that R tends to 


? The procedure applied in calculations determines magnitude of the above integrals but not the 
sign. In the formula for S we used the sign corresponding to negative [core i.e, to the lowering of the 
energy of the state caused by the inclusion of the nuclear charge effect. 
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zero and the charge inside the box is distributed over its cross-section with a constant 
density. We shall investigate all the excitations which have been discussed previously 
for the above molecules (see I — for ethylene and butadiene, Olszewski 1957a — for 
Simpson cyanine dyes — henceforth referred to as Il). Thus we restrict ourselves to 
the z-electron systems. In the calculations of the states we include, as before, 
the interaction between the electrons. In the case of ethylene and butadiene we 
consider the configurational interaction in our computations. 

The problem of the mutual interaction of electrons in the model has already been 
discussed earlier in detail (see I, III). The interaction integrals and energies corres- 
ponding to the present case are identical to those of the papers referred to above. 

It seemed useful for several excitations to present separately the relative displace- 
ment of the states caused by the nuclear charge effect with respect to the ground 
state. 

The question of the choice of Z for the calculations furnishes another problem. 
In the case of the z-electron approximation, e. g. for polyenes, we usually assume, 
as mentioned at the beginning, Z equal one. In this way we assume that the system is 
electrically neutral. However, the calculations of the nuclear charge effect for ethylene 
and butadiene already indicate that for such Z the difference in excitation energy as 
compared with the simple model or one including only the mutual interaction of the 
electrons is very large, and at the same time the excitations energies themselves 
become extremely high. ! 

Simultaneously, the configurational interaction begins to play a very important 
role. The depression of the ground state caused by fhis interaction raises the excitation 
energies still more rapidly. The configurational interaction effect ceases here to have 
a perturbational character, and the results become doubtful in the sense, that in view 
of the importance ot the effect, very many configurations should be included in the 
interaction. This, however, makes computations very difficult. For these reasons we 
must find a more reliable way of choosing Z, namely, we take the screening effect 
into consideration. ; 

When the influence of the nuclei on the eigenstate of a given electron is investiga- 
ted, the remaining n, — 1 electrons (x-electrons) screen the nuclei. This screening 
will be stronger in the case of a very thin stright tube, since, if we momentarily treat 
the electrons as point charges, almost the whole of the remaining electrons will appear 
on the line of interaction of a given electron with the nuclei. We assume approximately 
that the average distribution of these remaining electrons is such that it corresponds 
z . This corres- 
Te 
ponds to an equal accumulation in the centres of the Greene negative charge distri- 
buted along the entire free-electron path. At that time, for instance, for butadiene, 
instead of Z = 1, we take Z,g =}; for ethylene Ze = $. (One might also try to 
obtain the “true“ values of Z’s for the systems investigated, e.g. for polyenes, by a 


semiempirical method). 


n 
to a decrease in each of the charges of the n, positive centres by ~ 
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In the case of symmetrical cyanines (see Fig. 2) the choice of Z based on the 
principle of electrostatical balance will give: 


= 4(+) + 3 )-( = 1 (+) — for C atoms, 
Zy = 5 (+) + 3 (-) = 2 (+) —for both N atoms, 


where the removal of the z-electrons from the core and the presence of the (+) 
charge of the whole cation on the N atom is taken into account. On the one hand, 
it appears that it is not possible with our present form of the model (R — 0) to perform 
the entire computations when ہے‎ 4 Zc, on the other, we may assume that the (+) 
charge of the cation “wanders“ across the whole chain from one nitrogen atom to 
another. Thus we assume, approximately, that the whole positive charge becomes 


2 
uniformly distributed over all centres. Therfore Z = (2 a a + (n, — the number of 
c 
centres). As we have n, = n, +1 a-electrons in the cyanine chain, the positive charge 
Ze on the centres, after the screening is taken into consideration, is 1 + 
80 0 ۳۷ ms 2 
Me ne n, 
One might also assume that the (+) charge of the cation lies outside of the free- 
-electron bond path. Then, by an equal distribution of the positive charge over the 


1 1 
centres, Z should be (2 ا‎ , and Zæ, — —. This last assumption seems, howe- 
n 2 n 
c c : 
ver, to have a rather factitious support. 


Similarly, one might consider Z for the case in which the each cyanine has a loca- 
lized electron pair at the nitrogen atom (n, = ہہ‎ —1). Then 


Zy = 5 )+( +5 )-( = 0 — for one N, 
Zy = 5(+)+ 3(—) = 2(4+) — for the other. 


In view of the equivalence of both N atoms in the conjugated chain, we assume as 
the actual Zy the average of both values. We obtain the average Z, = l, — a quantity 
identical to Zo. It is characteristic, that just for this case we obtain for À the values 


i 2 
most close to the experimental data. (we then have بر‎ again equal to = ; 
c 

As in the case of n, = n, +l, we may assume also in the present case that the (+) 

charge responsible for the cation character of cyanine lies outside the free-electron 
path. Then Z = =, E PCR ای‎ 
Ne Ne Ne 

Up till now our calculations have involved constant electron density over the 

cross-section of the box (in the case of R tending to zero). It may be easily proved 

that if we perform calculations (for R > 0) with the actual radial functions of x and ø 


electrons taken according to (3), we obtain, owing to the normalization of these func- 
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tions, the same values for the differences of the [im integrals, or for the ت7‎ integrals, 
as in the case of the constant radial density. Yet, when the nn cn hice meen 
integrals are computed, the integrals calculated over the radial functions of the 
o-electrons alone (with the same condition on R as above), give the same values as 
when calculated by replacing the radial function by a constant; however, in the case 
of the “pure“ x — x interaction the interaction integrals will be 1.5 times greater 
than that calculated with F(r, 9) = const. This causes a corresponding change in the 
state system: the singulet-triplet distances, for instance, increase 1.5 times in compa- 
rison with those presented in (I). 

All the versions of calculation discussed above were investigated numerically 
for several molecular examples. The results are presented below. 

Ethylene. The energy of the lowest excited state is expressed by the formula: 


E (V,T) = F 1E ور‎ 6 [5 + Jus Se Kiz» 


where ,رط‎ E, are the eigenenergies of the states 1 and 2 in the zero-order approximation 


(see e. g. IT, eq. (1)), J, K — Coulomb & exchange integrals. The ground state has the 
energy: 


E (1) = 2E, +26 + Ji 


The relative change of energy of the (V, T) state with respect to the ground state 
caused by the inclusion of the nuclear charge effect amounts therefore to (in eV): 


ee کے‎ 5.80 Z 
For the state 2 (see I, both electrons are then raised on to the level of n = 2) this 
change is (in eV) 
"۲۷۶۷۹۹۶ = 11.59 Z 
Tables I, Ia present the diagram of states for ethylene. In this case there is no contri- 
bution from Z to the non-diagonal element of the configuration interaction matrix. 


Table I 


The electronic energy levels for ethylene in eV (without configurational interaction). 


AS FEMO method nuclear charge effect 


AS FEMO melded Exact LCAO 
method without = MO method 

ie nucl. charge „pure %—7| (Parr & Craw- 
effect interact., ford 1948) 


Zea = $ 
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Table Ia 


The electronic energy levels for ethylene (in eV) (configurational interaction included) 


AS FEMO method nuclear charge effect 


AS FEMO included Exact LCAO 
MO method 
State | method nucl. pure“ 7—7 Exper. 
h ffect 8 (Parr & Craw- 
charge eïtec Z= Zef = 4 interact. ford 1948) 
Left T 4 
0 0 0 0 0 
12.24 9,34 9.15 3.1 6.4 
14.10 11.20 11.58 11.5 7.6 
26.56 20.78 21.45 15.0 


Butadiene. The relative change of the energy of states caused by the nuclear 
charge effect (before the configurational interaction is included)is: (eV) 


for state Vas» Tog: 133 — Iz = 2.00 Z 
for state Viss T14: Jag — A = 8 Z 
for state Vas, Toy: JG — Bo = 3.327 
for state Vis, Via: IŞ — lir = 6.06 Z 


In the configurational interaction (of the one-electron excited states) in the matrix 
elements (Pol Vos) and (Ml V3) there appear the following expressions representing 
the nuclear charge effect (eV): 


The full set of elements are: 


(Vol Vis) = 6.55 eV, (Zola) = 10.22 eV for Z = 1; 
(Fol V8) — 2.23 eV, (Vol Vis) = "3:79 eV for Z ='4; 
(7 ol Vis) = 2.64 eV, (7 ol Via) = 462 eV for Z = ł, and pure 


7-7 interaction 
Tables II, Ila give the systems of butadiene energy states. 
Simpson’s symmetrical cyanine dyes. 
+ 
Me? N = CH — (CH = CH), — N Me? = 
+ 
= Me? N — CH = (CH — CH), = N Me? 
Fig. 2 


We consider here the locations of absorption maxima, which have been discussed 
previously for these molecules. They correspond to transitions from the highest 
filled to the lowest unfilled level. A good agreement with experimental data on thed’s 
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Table لا‎ 


The electronic energy levels for butadiene (in eV) (without configurational interaction) 


ee . 
AS FEMO method nuclear charge effect 
AS FEMO 


included Exact LCAO 
ss method without + MO method 
nucl. charge مس‎ 70-7 (Coulson & 
effect Za 1 Zefit = ł interac., Jacobs 1951) 
Left Ta + 
1A, (Vo) 0 | 0 0 0 0 
1B, (Vag) 535 735 ` 5.85 6.61 5.83 
` A, (Vis) 7.28 13.34 8.79 9.37 8.16 
1A, (Vox) 10.17 13.49 11.00. E 0 6.15 
1B, (Via) 13.44 20.82 15.28 16.25 11.51 
3B, (Tos) 3.88 5.88 4.38 4.41 —0.78 
SA, (Tis) 1223 13.29 8.74 9.29 2.76 
45 (وو1)‎ 9.94 13.26 10.77 11.14 4.66 
3B, (Tı) 13.31 20.69 15,15 16.06 7.26 
N —0.32 —7.27 —1.67 —2.28 —2.79 
Table Ila 


The electronic energy levels for butadiene (in eV) (configurational interaction included) 


AS FEMO method nuclear charge effect 
included 


AS FEMO Exact LCAO 


method MO method 
State without «pure nn Coulson Exper. 
nucl. charg. interact. Jacobs 
effect Za t 1951) 


corresponding to transitions mentioned above was obtained when only the interaction 
between the electrons was taken into account. Now, if we include the nuclear charge 
effect (with screening) we obtain in the case of n, = n, + 1 the values of A which 
differ in some points rather much from those obtained from the previous theory 
and from experiment. The present values show some displacement in the direction 
of the shorter waves(see Table III). On the other hand, in the case of n, = ہہ‎ — l we 
obtain a À somewhat, longer than that obtained from the previous theory and from 


a 
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Table III 


The positions of absorption maxima for Simpson symmetrical cyanine dyes. (ng = ne + 1) (in A) 


Amex Amax 
Mes AS FEMO AS FEMO method nuclear 

k n sais ‘ method charge effect included 

‘ model without nucl. 

charge effect 

0 4 2070 ` 
1 6 30 
2 8 45 
3 10 5880 


Table 1118 


The locations of absorption maxima for Simpson symmetrical cyanine dyes. (n, = n, — 1) (in A) 


ler 
Me, کت‎ AS FEMO method nuclear = 
ax ۱ 
7 simple method charge effect included (McGlynn 
model without nucl. & Simpson 
charge effect 1958) 


The A calculated here correspond to the excitation energies equal to the centre a gravity of the singulet 
and triplet excited states. They differ somewhat from the A calculated as the arithmetical mean of the 
singulet and triplet wavelengths, as was done before (see II) in the case of the FE states including only 
the electron-electron interaction effect. 


experiment. However, in this case the agreement with experiment is rather better 
than in the previous one. This is so especially for 4’s obtained if the presence of the 
(+) charge on the one of the nitrogen atoms is taken into account. (As it was men- 
tioned above, we then have Zy = Ze). The constancy of the vinylene shift is generally 
preserved in all the calculated results. 

In view of the above it seems that the localization of the electron pair on one of 
the nitrogen atoms of the chain is now established. 

It is possible in the case of symmetrical cyanines to calculate the absolute 
influence of the nucleus placed in the middle of the free-electron path on several 
states (with a paired quantum number n). For example, in molecule k = 0 (see Fig. 2) 
the energy of the state with n = 2 is reduced by 12.53eV (Z = 1), the energy of 
the state n = 4 — by 16.01 eV, of the state n=6 7 18.08 eV, etc. 
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It is very characteristic for calculations for ethylene and butadiene that results 
obtained with a screened nuclear charge effect gave for the lowest singulet excited 
states values almost identical to those obtained by Parr and Crawford (1948) and 
Coulson and Jacobs (1951) by means of the LCAO MO method. However the agreement 
with experiment is now, of course, worse than for results obtained in the free-electron 
model which included only the mutual electron interaction in the calculations. 

In general we conclude that the inclusion of the interaction effect between nuclei 
and electrons increases the distance between the free-electron states characterized 
by different quantum numbers n, and the greater the value of Z, the stronger the 
effect. 


§ 3. Concluding remarks 


The calculations of the influence of nuclear charges on the free-electron model 
eigenstates could be carried out for a three-dimensional box without any restrictions, 
i.e. we may imagine an arbitrary displacement of the nuclei inside a box and an arbi- 
trary choice of charges on these nuclei and perform calculations in each case. This. 
therefore suggests the possibility of investigating the way in which the optimal distri- 
bution of the nuclei depends on the magnitude of the nuclear charges, the number 
of electrons in the model, the box dimensions, etc. This situation, however, changes 
completely when we undertake calculations in the case of R — 0. Then the assumption 
of a constant distance between the nuclei along the free-electron path becomes essential. 
In earlier computations, of the kinetic energy and the interaction energy between the 
electrons, this assumption did not play an important role, besides establishing the 
total length of the free-electron path. Now, however, by any differenciation of atom- 
atom distances (or the atomic charge Z), the eqs. (12) and (14) are no longer fulfilled. 
It is then impossible to calculate the energies of the excitation states, In this situation 
it becomes necessary to return to calculation for a box with finite transversal dimen- 
sions. : ۱ 

It seems, however, to be of particular interest, that even if a constant distance 
between atoms is preserved (in the case of R tending to zero), not all excitations are 
possible. For instance, in butadiene (four centres problem), besides the excitations 
given in Tables II, Ila (1, 2-< 3, 4), the following ones are permissible: 


ko gana 8, 9 AT 122 13, 4% 


as well as all other transitions between the states of all these groups and in each of 
them. Yet transitions from every one of the states mentioned above to the states of 
quantum numbers n = 5, 10, 15 etc. cannot take place in view of the infinite transition 
energy. Similarly, in the case of ethylene, according to (12), the aliowable transitions 
between states are given by quantum numbers n grouped as follows: 1, 2; 4, 5; 1, 8; etc. 

Thus, in the system of the free-electron model eigenstatés, after the inclusion 
of the interaction between the electrons and the equidistant centres, a zone structure 
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appears. The number of states in each zone. is equal to the number of centres. 
The zones are separated by the forbidden states, in the sense discussed above, with 
quantum numbers equal to multiples of the number of the centres increased by unity. 

This effect is quite analogous to the Brillouin zone structure effect of the band 
theory of solids, which might have been expected. A more extensive discussion of the 
whole problem lies, however, beyond the purpose of the present work. 

lote added in proof. As mentioned above a very good agreement between 
calculated and experimental locations of absorption maxima may be obtained 
assuming n, = n, — l. The question of the number of z-electrons in the chain 
is connected with Kuhn’s rule for the effect of replacement of the middle CH group 
in the chain by heteroatom. 

It was namely clearly established by Kuhn (1951), basing on experimental 
data also, that when there is an even number of the double bonds in the limiting 
structure of symmetrical cyanine then the replacement of the middle CH group 
by the nitrogen atom causes the shift towards the shorter wave lengths; in case 
of the odd number of the bonds the maxima shift towards the longer waves. 

This rule may be qualitatively easy verified by our presènt scheme, too. 
Replacing in the middle centre Zo by Zyn, and taking into account that the 
nitrogen atom is more electronegative than the carbon one (Zy > Zc), the exci- 
tation energies tend in the case of the odd number of double bonds towards the 
higher values (the hypsochromic shift) and for the even number of double bonds 
in the chain — towards th lower values (bathochromic shift). This is in fair 
agreement with Kuhn’s rule, Such result may be obtained, however, in our scheme 
` only assuming n, = n, + 1. Thus it seems that in further investigation the last 
relation betwen n, and n, should rather be accepted and not the previous one: 
n, = re — l. 
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The Coulomb and exchange operators occuring in the Hartree-Fock equations for 
molecular orbitals of the linear free-electron model are given for the case of a potential 
box radius tending to zero in the form of explicit and finite functions of the position 
of the electron on which the operators actually act. The Hartree-Fock differential-integral 
equations are then replaced by ordinary differential equations. 

This has been done with the reservations that the Coulomb and exchange operators 
were calculated by assuming a suitably modified expression for the mutual interaction 
between electrons and an exchange potential of the Slater type averaged over the entire 
region of the potential box was introduced as the exchange operator. The calculation gave 
expressions reproducing the energies of the respective states, whereby the approximations 
introduced are taken into account, plus quantities that are the same for all states. The 
form proposed for the Coulomb and exchange operators.is that obtained from the above 
calculations, where those parts contributing to the reproduction of the above-mentioned 
uniform quantities which exceed the energies of the calculated states, or which gave 
a value of zero, are discarded. | 


§ I. Definitions and basic formulae 


If the solutions of the Schroedinger equations for a given one-electron problem 
are known, and if the electrons to which these equations apply form a multi-electron 
system, then we can write equations whose solutions will be better in comparision with 
the previous eigenfunctions and eigenenergies of each electron in this system. These 
are the Hartree-Fock equations which take into account the interactions between the 
electrons and the antisymmetrization of the wave function of the entire system in 
accordance with the Pauli principle. 

Mathematically, the Hartree- Fock equations are integro-differential equations of 
a rather complicated form. It i is especially not clear how the independent variable of 
the function being actually sought is involved in the Coulomb and exchange operators 
entering into the 7: 

(121) 
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There arises the question whether it would not be possible in the case of calcula- 
tions for one-electron orbitals which are solutions of the eigenproblem of the free 
electron model and for the limiting case of a very thin potential tube to give the inter- 
action operators of the Hartree-Fock equations in an explicit and finite form. At the 
same time the integro-differential equations would be replaced by ordinary differential 
equations. 

It turns out that this is in fact possible, however with certain restrictions. The 
aim of the present paper is to show the way in which the required form could be 
given for interaction operators in the above-mentioned case. 

The Coulomb and exchange operators, as defined in the case of molecular 


orbitals by Mulliken (1949) and Roothaan (1951) have the form 


gi) =e ( f: oe) ar] pi) (la) 


|r 
: r) g(r 
وہ‎ = à ( f EES ar) وم‎ (1b) 


In addition to the conditions of linearity and hermicity, these operators should, 
of course, satisfy the basic relation 


fg (Dg; oy (0 dt = Jy = | o (© gj; () dr (2a) 
[GO Ro) dr = Kg = | G0) Kr) dr (2b) 


where J;;, K;; have the conventional meaning of Coulomb and exchange integrals. 

This definition obviously could have been extended automatically for the case 
of our FE model by putting for g functions that are solutions of the Schroedinger 
equation for a particle moving in a cylindrical box and which satisfy the proper boun- | 
dary conditions (see e. g. the previous work of the author, hereafter denoted as IV). 
This seems, however, to be of little practical value because a considerable compu- 
tational difficulties. 

We define, somewhat similarly as above, the Coulomb and exchange operators 
for the case of a box whose cross section tends to zero, where this limiting assumption 
is already incorporated into this definition. Thus, in case the cross-section of the box is _ 
very small we can replace, by way of approximation, the interaction between charges 
whose radius vectors have projections on the z axes z and 2’ by the interaction of two 
charged cross sections of the box at positions z and z’, perpendicular to the axis, 
where the charge is assumed to be uniformly distributed over the entire cross- 
-section. 


The operator of the mutual interaction between electrons in positions T (z,r,@), 


1 1 
r’ (z’,r’ 0’) is then defined, instead of mr] as the averaged r-r] over the whole 
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interacting cross-sections: 


R R OR OR 
2 
E Zy z = E á 1 f / 
(z, 2’) (RI? JIT fr pew dû dû (3) 
0 4)4 ٤ 


R — the radius of the cross-section of the model box. 


Let us now define new Coulomb operators 4° and exchange operators K’ by means 
of the operator E (z, z’)!: 


Ga Y z) = lim 5 ZA Z z’ 7 z z 4 Aa 
1 (z) [ m [GE €, ) y; (2) d [۷ (ھ)‎ ( ) 
oe yY Z) = lim ا‎ z E Z z y A dz’ A Z 

1 ( ) [ zi i 1 ) ( 2 ) ( ) ] 1 ( ) (4b) 


where y is the part of y dependent on z. We define the Coulomb and exchange integrals 
derived from the operators (4) 


j vî ei vw (de = J; (5a) 
[¥ORxOu=K; (5b) 
: a à : A eee 
Operators (4a, b) are linear and hermitian, since the functions y, = [| “72 


are real (L — the length of the model box). From this and from the definitions (4a, b) 
it follows that for 2 an equality analogous to the equality on the right-hand side 
of Eq. (2b) will also be satisfied. We shall show that also in the case of J;; an equality 
corresponding to the right-hand side of equality (2a) is satisfied. But the integrals 
ee K;; are not identical to J;;, K;; for R — 0, since the procedure of passing to the 
limit in the first case involves the integral occurring in the definition of the operator 
and in (5a, b) under the integral sign. 


The entire scheme presented above will be useful if the defined operators can 
reproduce well the actual values of the Coulomb and exchange integrals. In other 


words, it should be investigated whether and by how much the: valeus K;;, J; differ 


from K;;, J;;. To do this we must calculate the integral in (4) and then in (5). The 


integration in (4) gives us g’, K’ as an explicit function of z. 


As has been done earlier (III, IV), we replace in E the operator by an 


ly — r'| 
1 

integral representation of the Green function G of the three-dimensional Laplace- 
-Poisson equation (see Olszewski 1957, referred to hereafter as III, Eq. (12)). Also, at 


E SS ee‏ ھچ رہ عتت 


A similar operator is given by Ham and Ruedenberg in their semi-empirical method of treating 
the electron interaction in the FE model (1956). 
# 
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the same time we transform the product y’;;(z’) رپ‎ (2) so that each integral of the 
operators g’ and &’ in (4) can be split into two integrals of the type 


e2 


= lim f+ cos cz’ dz’ E (z, z') dz' = lim — l p | co یں‎ An G (r,r) do do’ dz’ 
R-0 مھ‎ L (aR)? 

(6) 
where do = rdrd#, and the integration over dz’ extends from 0 to L (we integrate 
first over z' and then over the variables of integration of G, what is permissible, since 
as in IV (Olszewski 1959) we can show that the integral over x, k is uniformly conver- 
gent for the entire range of variability of z). As a result we have (case cA 0) 


C= 7 ‫َ {cos cz Ci (cz) + cos c (L — z) Ci [c (L — z)] + sin cz Si (cz) + 


+ sin c (L — z) Si [c (L — 2)]} + @ (7) 
wheré 
2 
a= — Ê «20005 | - lim log R + > + log 2 — y — loge | (8) 
= R-0 


y is the Euler constant. Besides @! we shall also calculate the operator @, corresponding 
to the limit lim @: 


0ء 


2 : 
es = — S logs + log E — 2) +2 ( — im log R + À +1062) (9) 
We can now calculate the mean values of the operators Fy ; K; (through @, and یھ‎ 
according to (5): wé multiply by suitable y and integrate over z. We obtain 
a combination of integrals of the type (cf. (16) and (20) in I — Olszewski 1955)): 


fe cos az dz; fe cos azdz; Je dz; fe dz; (10) 
unfortunately, it turns out that the value of the integral 
= Je: cos 52 (11) 


coincides with C obtained in I, eq. (18) (see also Olszewski § 4, 1957) only if c Æ a.? 
If this condition is not satisfied C’ is essentially different from C. Thus it is seen that 


2 In a relatively simpler form than that obtained in I (after applying the boundary conditions of - 
the linear FE model) C is given by 


ir tiga g 
= cos cz’ cos az : dz dz = 
Rao rR یت‎ 


1 : T 
= EE ES QE my} ا ای‎ Sie FM S EMI} انا‎ 


where cand a are n+m, n—m, from Jpm» Kpm respectively. If the condition c a is satisfied then we at 
once obtain from the integration in (11) the result (æ), where C is already supplied by the first term in 
expression, (7) for, Êi Q: gives zero because of the orthogonality of cos cz + cos az. | 
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the interaction operator we have defined is not in a state to reproduce the exchange 
integrals and integrals of the type J;;3 we must seek other methods to overcome this 
difficulty, e. g. approximation methods. On the other hand 


| az @ dz = Si (aL) - te 
a 
correctly reproduces the corresponding element of the Coulomb integral. The same 


C 0 2 
applies to f @,dz. We must still investigate how the operator @, reproduces £. To do 


this we write the operator somewhat differently than in (9): 


; e? 
@n = — + [log z + log (L — 2) — 2logL — 1] + O (9a) 
where 
; ice : 3 
Qo = — 2— | log L — lim log R + — + log2 (9b) 
L R+0 4 


(à ۰ f ۰ ۰ 1 7 
The part of the operator وھ‎ without @, gives, after multiplication by 7 and inte- 


2 
3 à € > : : : 
gration over z, precisely T’ i. e. the value of interaction integral C when c and a 


tend to zero simultaneously. Thus, by summing all the respective terms fromm (10) 


we finally obtain: 


J = f YD Fi yed 
یں ما می 2( ۷ئ‎ tate 12) (12) 
R= 0 2 
= f| یی‎ (2) dz 


[2\f. na t 
for y, = (7) sin — 2. 
Hence our interaction operator reproduces the integral J;;; however an additional 
term independent of ij occurs here. Nevertheless, as will be shown below, this seems, 


in general to be sufficient for our purposes. 


It should be noted that our procedure (with the aid of the expansion of the Green function) of 
calculating the interaction operators (of the type  & @)) and previously of the interaction integrals 
(see I, III) fixes the value, but not the sign, of these quantities. In calculating the elements C of these 
integrals (see e.g. (x) footnote 2 of this paper) and the elements of the operators we put for these quantities 
a sign that corresponds to a positive sign of the interaction integral composed of these elementi or obtained 
from these operators. As regards the operators, this choice is unique for all @’, while in the case of the 
interaction integrals C with the sign given in («) (and in the previous papers, i. e. positive C) gives 
the plus sign for all Coulomb integrals and almost all exchange integrals. ia the case of the latter, an 
interesting exception is Ky, 9,41 When n > 2. Their absolute values are, in general, very small frac- 


tions of the integrals Ky. J 
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2. Approximation of the exchange operator 


The point of departure for the approximation of the problem will be the Slater 
modification of the Hartree-Fock equations. In this case the equation for (م) رہپ‎ takes 
the form: 


N/2 f e? 
H g; (r) = 1 2 Pr (T) pr (T) — Cav, ex (T, 27 ror] رپ‎ (r) = E;qi(r) (13) 
hol 


N is the total number of electrons, H is the one-electron operator of the eigenvalue 
problem. 


We assume here, in accordance with Slater’s treatment, that in each state k 
there are two electrons whose spins are oppositely directed. ہے ہہو‎ (1, T) represents 
here the so called exchange density averaged for all states k. The entire expression 
under the integral sign now has the standard form of an interaction potential of the 
Coulomb type, which corresponds to the action on the electron in position r of the 
total charge density of all electrons in position +” and corrected by the averaged ex- 
change density corresponding to a loss of charge of the same spin as the electron 
in question in its direct vicinity. 

In the case of a totally free-electron gas 0,,,, (T, T”) was first calculated by 
Wigner and Seitz (1933). It was then assumed that the electron momenta are directed 
in all directions. We shall now calculate Oy ہے‎ (T, T) for the case of an electron gas 
enclosed in a very thin tube of length L. We can then assume that only the momentum 
components parallel to the box axis are of importance. 


The exchange energy is 


1 aa / 1 7 
af [hr tol’) |? -|)۸| [*|(٭|۔ء‎ dr dr (14) 
where, if the electron densities for both spins are equal, 
N/2 š 
relr) = (rjo) = 2 pe (r) وص‎ )۳۰( - (5) 


For the functions y we now take the quasi-classical WKB approximations. Taking 
into account the fact that R (the box radius) tends to zero and that the electron momen- 


tum in the tube can be directed positively or negatively with respect to the box axis, 
we obtain in place of (15) ۱ 


P 
| | 
=P nu 
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(see Brillouin 1934, chapt. 7). We replace the summation over the states by an integral. 
P is the momentum of the electron in the highest occupied state equal to AN Final- 


ly, in the case of a very small R the expression for A takes the approximate form 


2 


7 sin > |z— 2’ | 
Ara ۹ی00‎ £ 2 h f 
72 am 3 eae u P dt dt (16) 
اوھ‎ 


where و‎ is the classical density expressing the number of electrons per unit length of box 


We regard the expression 
P 2 
sin — |z — 2 
f h 
Oav ex (z, ا‎ > 5 P É ; (18) 
7 rai اس‎ | 


as the averaged exchange density. Eq. (18) represents the electron density function in 
the hole formed in the vicinity of the electron at the position z in a linear distribution 
of charge of N electrons. In the direct vicinity of the electron under consideration 
(z’ = 2) the loss of the density is equal to half the total electron density, which accounts 
for the fact that another electron of the same spin cannot be in the neighbourhood 
of a given electron. As |z—z’| increases, 0, (Z) 2’) gradually decreases until 
a charge density equal to the total electron density is attained. 

According to Slater, in order for the averaged exchange potential to cancel out 
the consequences of the fact that in the Coulomb potential the terms corresponding 
to the self-interaction of the electron was taken into account it is necessary that the 
charge in the exchange hole equal the charge of the electron. Indeed, in our case we can 
simply represent the magnitude of the charge in the hole as a function of the maximum 
momentum P and the actual position on the tube axis. It is thus seen at once that for 
a large number of electrons (large P) and for positions not too close to the ends of the 
tube the charge in the hole is le. (In the case of the centre of the hole being very 
close to the end of the tube, this quantity approaches 1/2e). 

Except for the exchange integrals, the operators of $ 1 still do not reproduce 
the integrals J; = K;,; we must therefore also approximate in Eq. (13) g} (rp) 
for the case k = i. Here we set the electron density averaged over the entire area of 
(classical density). Thus, instead of 


1 
the box equal to IRL 


N/2 A | 
2 D où (^) یم‎ (r') (19) 
k=1 
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we take 
1 
2 A Daor 2 pare (19a) 


In order to decrease the error caused by the last approximation (see also the note 
below the Table I) we discard in the Coulomb potential in (13) the term corresponding 


Este 
x R?L 


It is necessary to take into account this correction in the expression for the averaged 


to the self-interaction of the electron. Then in (19a) it should remain only one 


exchange density and exchange energy. That is why in the product |(r| o, |’) |? in (14) 
we replace one of the (r|o,|17’) by 3 


N/2 
) 1 
mjo |y =A Pr (1) pr (r°) ا‎ (20) 


For the exchange density we then obtain‘ 
sin 2 |z—2'| 
h 


P 


e‏ ا 


SRP ; 
sin = |z— z'| 
IF , 0 h 
Oau. ex. (20 = 9 


(21) 


=| 


We now make a further simplification of the calculations with ọ'(z, z’) by averaging 
this expression over the entire region of the free path of the electron: we integrate over 
the entire region of the exchange hole’at the position z, then change the position of 
the hole from 0 to L. We obtain 


1 =; ; 4 1 
7 Wa = Cowie (z, 2’) a یڈہ‎ À 108 (Nn) + Ci (Na) + 


1 : AN N 2 
+ Na Si (Na) — nsi 0 3 — y — (cos À x - i) il l (22) 


The subscript of W indicates that the averaging refers to an object whose states 
from 1 to N/2 are occupied. | 


* It would be correct to neglect the term g% (r’) ہم‎ (7’) in 2 3 pr (r’) g (r) OF (r) Pp î instead of 
ee) 
cancelling yj (r) p; (r) by means of aa in (20). The error in the case of the FE function does 
not appear to be essential, especially since we take further the average of Qav.ex Over the entire range 
of both variables. ۱ 
* Despite the fact that we neglected the self-interaction of the electron, the volume of the exchange 
hole decreases only by 1/Ne (for large N). This follows from the fact that we obtain the expresssion for 
Cav. ex. (2,2’) by averaging over all states k. 


On the Coulomb and Exchange Operators 129 


Thus for small R Slater’s equation now takes the form 


Hy; (r DRT AE 76+ ' o; )۳( = : 
gi (r) + À FE) mr) RL (rm) = O) 
Ri 


where we have omitted the term k = i in the summation. 

From the point of view of the suitability of the approximations (19)—(22) it is 
essential to check how the exchange energy and the energy of the total electron inter- 
action calculated by means of these approximations differ from the corresponding 
energies in the Hartree-Fock scheme. Table I shows us these quantities calculated 
in both ways for the ground states of several multi-electron systems. 


Table I 


2 
e 
Electron interactionen energy in units of I for systems of N electrons 


Exch. energy Total interaction 


Total interaction 


orna a Hatiro es energy of electrons | energy according 
by Ai exch. energy approx. to Hartree-Fock 
0.366 0525 RE 3.54 237 
1.255 1.265 8.50 7.22 
2.178 2.266 16.8 15.2 
3.124 3.434 28.4 26.5 
12.88 21.41 350 339 


The approximated exchange energy is calculated in Table I by taking Eq. (16) 

with 0... (2 2’) instead of @,,, .,(2, 2’). Similarly, for the Coulomb energy of the system 

2 

derived from the approximated terms (J;;) we can set a 7 . The exchange energy 
calculated according to the Hartree-Fock scheme given in the table does not of course 
contain integrals of the type K;;, because of the discarding of the term k = t in the 
exchange operator. It is characteristic that if the total exchange energy including K;; 
is calculated, and therefore with the help of @,,, ex (Z: 2’), the value then obtained is 
too small to take into account the cancelling of the corresponding terms of the Coulomb 
energy by the exchange energy for the case k = i. 

The results of the calculations presented in the table show that the error in | the 
calculated exchange energy increases considerably with an increase in N, while the 
error in the total electron interaction energy distinctly decreases with an increase in N, 
for then the Coulomb energy plays an ever more important role. We can treat this 
result as a mean result holding for all states i occupied by N electrons. 
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Comparison of the results of the total interaction energy with the accurately 
calculated results (in the Hartree-Fock scheme) suggests that for large N we would 
obtain a better approximation by simply setting W y equal to unity. In the case of 


2 
very large N, the value Wy indeed tends to 1. Then Eqs. (23)—(26) undergo further 
2 

simplification. | 

The use of equations (23) assumes practically the same exchange energy for all 
states i of the system. This assumption is justified since the Hartree-Fock exchange 
energy has relatively small influence on the energy of the states of the FE model 
(and also for Ai not too small, on their difference), and the more electrons in the 
given system, the less the influence. 


§ 3. Summary 


Let us now return to the interaction operators defined in § 1. If we replace the 
operators occuring in (23) and therefore of type (1), by new operators (4):. 


2 
Fv; ھ)‎ = Hv; (2) + { È پ27‎ + (1 — Wn) Cp} ¥; (2) = Ej; (2). (24) 
bay i 
and with their help calculate the energies of the individual states E; with the functions 
w; (2) in the zero approximation, i.e. with the FE MO functions, then it turns out that 
for each i we obtain E; equal to E; calculated from equation (23) for the case R — 0, 
plus a quantity which is the same for all states. 

In the majority of multi-electron problems, however, we are generally interested 
not in the absolute position of a given state, but in the relation between the states of 
the system, especially in the correlations between the states. (These correlations are 
given by the Coulomb and exchange operators). The constant E occuring in 
all states do not then. play any role. 

On the other hand we can make for all states a unique separation of each of the 
operators J’, @, occuring in Æ so that one part of these separated operators inserted 
into (Æ will give, in the case of the zero approximation functions y; (z), the energies 
E;= f y; (2) F y; (z2) dz identical to E; in Eq. (23) (R > 0), i.e. the energies of the states 
with the effect of the remaining electrons taken into account. Hence, this part of the 
operators describes the correlation of a given state i with respect to the remaining 
states. On the other hand, the remaining part of the separated operators (obtained 
from $’, 0, being of the same form for all states, when inserted together with the 
first part into (Æ gives the energies E; = E, +B, where B is the same for all states. 

Because of the correlations interesting us, it therefore seems proper to discard 
in (24) the components of the operators f’ and 6 which contribute to the production 
of quantities that are the same for all states and which exceed the energy of these 
states (shifting the entire spectrum of these states), or which give a contribution to 
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the energy of the state equal to zero. These will be O; and O respectively. In other 
words, it should be expected that after discarding the above mentioned parts of the 
operators ء2‎ (in Q’) and @, as capable only of shifting the system of states as a whole, 
without affecting the mutual positions of these states, the obtained interaction operators 
properly represent the effect of the remaining electrons on.the localization of the given 
electron state among states corresponding to the remaining electrons. | 

Inserting for J’ the differences i -م6)‎ €) and neglecting the above-mentioned 
terms, we obtain for Eq. (24) the form 


N/2 


2 
H y; (z) + D (— 1) i [log z + log (L—z) — 2 log L— 1 — cos cz — 
k=1 
ki 


— cos c(L— z) Ci [c (L — z)] — sin cz Si (cz) — sin c (L — z) : Si [c (L— 2)] | + (25) 
2 
OM z + log (L—2)— Plog L— 1} wo) = 
2 


= E;y; (2) 


where 
= — (26) 


In this way we have replaced the integro-differential equations of Hartree-Fuck 
for the case of the free- electron model, by means of the modification of Slater, by 
ordinary differential equationsÿ. 

The operators corresponding to the Coulomb and the averaged exchange 
operators can be represented here in explicit and finite form, valid, in general, for 
all linear systems of electrons with even spinsf. 


5 If we were to define the Coulomb and exchange operators only as operators satisfying the properties 
of hermicity, linearity, and able to reproduce the Coulomb and exchange integrals, then in this case we 
could “guess at“ exchange operators with these properties. They will be the differences of the type @} 
with OL discarded. With such a definition we refer to the relations between the operators J, K and the 


mutual interaction operator (through E(z, z’), if only these relations help us guess the form of Z; 


r— r'| 

and K;. Such á procedure might be ambiguous, however, and that is why we keep the operators Z’, © 
calculated by means of Eq (7)—(9). Finally, as the interaction operators we introduce such 3” as differ 
from (f’ only in the discarding of terms giving the previously mentioned factors constant for all states 
and exceeding the energy of these states. 

6 We have not obtained an here operator with properties analogous to F representing the influence 
` of the nucleus on the eigenstates of the model (H, as an operator of the Schroedinger equation of the one- 
-electron problem, is only a kinetic energy operator). In the case of many electrons it can be assumed that 
because of screening (see Olszewski 1959) the effect of ص2‎ (effective nuclear charge) on the relative 
position of the states is slight; we again have here, in a rough approximation, a 6 of the entire 
system of states without any essential changes in their relative positions, Only the division of states into 
zones owing to the presence of positive centres seems to be especially significant. 
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The form of equations (25) may prove to be useful especially for a simplified 
solution of the Hartree-Fock equations in the case of the multielectron problem in 


the linear FE model. 
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L'INFLUENCE DES MOMENTS QUADRUPOLAIRES SUR LA 
CONSTANTE DIELECTRIQUE D’UN LIQUIDE DIPOLAIRE 


par 
A. R. 7 
Institut de Physique Théorique de L'Université Adam Mickiewicz, Poznan 
(Manuscrit regu le 24 juillet, 1958) 


La théorie de Onsager est appliquée aux liquides polaires dont les molécules présentent 
des moments dipolaires et quadrupolaires. Un modèle ayant la forme d’un dipéle excentri- 
que dans une sphére entourée du continuum diélectrique est employé. Le calcul donne 
pour la polarisation un terme mixte supplémentaire dépendant des moments dipolaires 
et quadrupolaires. On montre que les valeurs des moments dipolaires pour certains com- 
posés chimiques liquides peuvent étre améliorées si les moments quadrupolaires ont 
été mesurés. Š 


1. Introduction 


Les travaux (Zwanzig 1956, Buckingham et Pople 1955a, Buckingham et Pople 
1957) sur l’influence des moments quadrupolaires des particules sur les propriétés 
diélectriques ne prennent pas en considération l’orientation des moments quadru- 
polaires sous l’action d’un champ électrique. Évidemment, le champ uniforme n’oriente 
pas un moment quadrupolaire. Cependant, le champ local, dont l’action sur une parti- 
cule nous intéresse ici, n’est plus uniforme dans les milieux denses par suite de l’action 
des particules proches. | 

Dans les travaux de Buckingham et Pople (1955b, c) les interactions moléculaires 
sont pris en considération. Une part de leur théorie concerne le gaz imparfait des 
molécules ayant des moments dipolaires et quadrupolaires permanents. Buckingham 
et Pople se bornent à l’étude du couplage deux à deux, c’est à dire à l’interaction de 
deux molécules les plus rapprochées. Les auteurs évaluent qu’il faut une valeur du 
moment quadrupolaire dépassant 20.10-?6 u.e.s. pour obtenir des résultats mesur- 
ables. La cause en est probablement la restriction que nous venons de preciser ci-dessus 
(d’ail leurs entièrement suffisante dans le cas du gaz), et, probablement, le fait 
que la polarisabilité des particules n’a pas été prise ert considera ion. Nous verrons 
plus loin que, dans le liquide, l’effet quadrupolaire est plus important que dans un gaz. 


(133) 
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On peut attaquer le problème de l’orientation des particules d’un autre point 
de vue. Considérons une molécule polarisable ayant un moment dipolaire et un 
moment quadrupolaire permanents. Le dipôle s’oriente dans le champ local; le moment 
quadrupolaire s’oriente aussi, ce champ n’étant pas uniforme. L’action sur le quadru- 
pôle est moindre, car le moment quadrupolaire est très petit; néanmoins ce dernier 
ne peut pas être négligé, vu que le Gradient du champ local est considérable. Nous 
voudrions démontrer que les composantes moyennes de ce Gradient ne sont pas 
nulles. 

La question de savoir dans quelle mesure l’orientation des molécules ayant des 
moments quadrupolaires est importante pour la théorie de la constante diélectrique 
forme le thème du travail présent. Vu les difficultés des théories statistiques des 
liquides dipolaires, le calcul se base sur la théorie de Onsager (1936) avec un 
modèle simple de la molécule. Il paraît que la méthode de Onsager appliquée au liquide 
pourra être utile malgré son manque d’éxactitude. 


2. Modèle modifié d’une particule de Bell-Béttcher et calcul de l’énergie totale 


Suivant Böttcher (1952, p. 135), nous admettrons qu’une particule est remplacée 
par une cavité sphérique de rayon a contenant un dipôle ponctuel m placé hors du 
centre de la sphère (Fig. 1). Le dipôle est polarisable; il ne doit pas être nécessairement 
dirigé suivant la direction du rayon de la cavité. L’entourage est un milieu continu, 
e étant sa constante diélectrique. 

Le potentiel ®, d’un ensemble arbitraire de charges ponctuelles situées dans 


une telle cavité a été calculé par Scholte dans sa Thèse de Leiden (citée d’après Böttcher 
1952, $ 16): 


co 1 
S 2 LE LM تی‎ pen) 0 
FR E 1-لسر ار‎ PO (cos 9) e”? , (1) 


P™ (cos 0) étant les fonctions de Legendre. C’est le potentiel dans la cavité dans un 
point ayant les coordonnées sphériques r, ®, ھ‎ Le potentiel du champ de réaction est 


coo 1 
D' = -D = - » 2 Cy سار‎ r PRE (cos 0) e’, (2) 
1=0m=— 1 
où . 
_(e-1) 0+) ہو‎ 
مہ‎ aa “4 
et 


1 
a bm, 7 PE: (cos #) ge (4) 


~ 
Il 


0 m= 


est le potentiel de ce même ensemble de charges placé dans le vide. 


L'influence des moments quadrupolaires 


3 s 
۹ ENS 


= se 


Fig. 1. Modèle d’une particule s — distance du centre au dipôle. 
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Fig. 2. Les axes des coordonnées x’, y’, z’ sont dirigées suivant les axes principales = LS ce Re 
iê laxe Z — suivant lardirection du champ extérieur, l’axe z — vers le point où est situé le dipôle. 
? 
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Pour le dipôle ponctuel déjà décrit, les b™, ont la forme 


bD; = z (m. + = m) si-1 (5) 


les autres bm, étant nulles. m,, my, m, sont les composantes du moment dipolaire m 
suivant les directions des axes des coordonnées x, y, z (Fig. 2), m étant situé dans un 
point (0, 0, s). Dans ce cas 


ge = - ¥ ر‎ r7 s1 [m, P® (cos 0) cos p + m, P®, (cos 0) sin g + lm, P% (cos 0]. 
i=0 
(6) 


En présence d’un champ extérieur E uniforme, la molécule est soumise à l’action du 
champ local F composé d’un champ de cavité vide E, 


3e 


hack 2e+1 pe 0 
et d’un champ de réaction R = — grad ©’, et on obtient de (4) et (5) 
R, = 3 لو‎ (1—72)# P® (t) cos g (m, cos p + m, sin p) — 
1=0 


(1) 
— (1—7)#7 ae cos @ (mz cos پ‎ + my sin œ) + 


1 . . 
ay (l= 72% PP, (T) sin p (my sin y — my cos p) + 


"02 
+ 12(1— 72)% PO, (x) cos p (2ج --1)] ے71‎ 7 — cos | 3 


Ke 5 ا‎ SE |: (1 — 13⁄4 P®, (x) sin P (mz cos p + my sin g) — 
1=0 


dp» 
7)0 ar a sin p (my, cos p + my sin gy) + 


1 $ , 
T (= 7% ام 3ء‎ (a) cos و‎ (my sin p — my cos p) + 


; ; . (0) 
+2 (1— ری‎ P (7) sin 77٦0 تع‎ sin y 7 à 
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R: = y Gi pik gE [i P™®, (7) (my cos p +m, sin g) + 


1=0 


ap 
۴)1 — 7°) oe gi (my cosp + my sing) + 


(0) 
+ 9م ہو‎ (Dm, +10 — 1) — m | A (8) 


où t = cos Û. Le champ de réaction d’un dipôle excentrique a été aussi calculé par 


Scaife (1953) et Frood et Dekker (1952). 


Le champ de réaction au point où est situé le dipôle a les composantes 


R, = fi Mys R, = famy, R, = fam; (9) 
où 
oo 1 co 
a a 9-لاو ری 2 کے‎ a ee 7 29-لاور‎ 12, (10) 
1-0 . 0 
Quand s = 0, c’est à dire quand le dipôle est central 
2)6 — 1 
h=h=h= چ‎ - (11) 


Quand s/az 0, un facteur de correction peut être introduit dans (11) pour obtenir 
la valeur de f exigée par les équations (10); quelques valeurs de ce facteur ont été 
calculées par Böttcher (1952, p. 138) pour f, (notation introduite dans le travail présent) 
et par l’auteur pour fı = fg (Tabl. 1). 


Tableau 1. 


Facteur de correction de fs 
e=2 


Facteur de correction 


| de fı = fz 


0,1 1,04 1,04 
0,25 1,05 
0, 50 S 1 ,67 


Un champ local agissant sur le dipôle potarisable, celui-ci acquiert la valeur m - 
dont les composantes sont: | 
| M = ka F Bay Fe == fj + Bay E, ا و‎ Ry,» À Sa if 2, 3, (12) 
où u représente 16 moment dipolaire permanent et a, le tenseur des polarisabilités, 


æ, F, signifie Dit a y = 1,2;3 etc. De cette manière le champ de réaction R, 
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à l'endroit où est situé le dipôle polarisable a les composantes: 
Ry =f; (u ap Bay Ry je Bry Eg)» A= 1, 2, 3. (13) 


Les indices À = 1, 2, 3 marquent respectivement les composantes x,y,z du 
vecteur. 

La part Wp de l’énergie totale d’une particule, liée au moment dipolaire (cf. 
Böttcher 1952, §§ 25 et 27), se compose de l’énergie W, de polarisation du dipôle, 
de l’énergie de polarisation W, du continuum hors de la cavité, et de l’énergie Wa 
d'orientation du moment dipolaire: 


1 


A=1,2,3 et » = 1, 2,3 etc. 


1 
W P a Patil, 


où Fa,,F, signifie 2 س‎ 7٦ 


yp? 


Une autre part de l’énergie totale est celle (Wg) d'orientation du moment qua- 
drupolaire dans un champ local non-uniforme F. D’après Rosenfeld (1951) 


Wo = - Qu (Grad F),, (15) 


Q étant le tenseur du moment quadrupolaire 


OET (16) 


où r détermine la position d’un élément de la charge électrique dans la molécule, et < > 
signifie une moyenne sur la distribution des charges électriques dans la molécule. 


Grad F est le tenseur dont les composantes sont (Grad F), = OÙ aa = %,; 


Fi 
ox, 
Xa = Y, X3 = z. Remarquons que Grad F = Grad R, E, étant uniforme. 

Avec les coordonnées x, y, z (Fig. 2) et la formule (6) nous pouvons calculer 


les composantes de Grad R. Le résultat étant bien compliqué, nous n’employerons 
que la moyenne d’espace, Grad R, définie par la relation 
a n 27 
Le 3 ۰ 
Av = Tr ۳ r? dr sin 0 dô dy. (17) 
0 0 0 


Dans les coordonnées x,y,z de la Fig. 2 les composantes moyennes du Grad R 
calculées avec Ø’ de (6) en tenant compte de l’orthogonalité des P™, (cos #) sont 


OR, OR, 14 dR, _ ƏR, OR, 

M یمور‎ bs plies ماس‎ Mae cuil 
OR, OR, OR, OR, 42 3 Ÿ 
Be de ap past (Sant x 


al—2 
ا‎ à (18) 
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Le procédé d’intégration est semblable à celui de Böttcher (1952, $ 16). Pour une 
molécule présentant les polarisabilités principales وہ ,وہ رھ‎ suivant les axes x, y, z 
nous obtenons d’après les équations (13) une expression de la forme 


ecu eal} 
Ra = ere n e ھ۶٦٣‎ (19) 
Alors énergie totale de la particule s’écrit d’après (14), (15), (18) et (19) 
ا ا ات‎ ee E as و‎ Qua | IRs 
Er از‎ i TS ےہ ور-1‎ 2% ur 


Hz j Qos «a OR, 
+h (ee ا‎ Sa Ga me) + 


Us sh 14 (Qui + Qo) a Qs3 %3 | 
Th (5 5 (1 — fy ag) ATEN TE as) 7 


1 a œ 0 
ESTA Wipe 7 و و جا 2 0ي‎ 
ous ) ۳ لفاون ےر‎ ree, ae) گے‎ +۶ 
+ termes indépendants de E,, (20) 


où ما‎ lg, وا‎ représentent les cosinus des angles entre laxe Z (direction du champ 
extérieur) et les axes principales du tenseur des polarisabilités. On peut choisir les axes 
x,y de telle manière que 4 = 0, J, = sin O, lg = cos O, où O est angle ZOz (Fig. 2). 


3. Calcul de la polarisation 


La valeur moyenne statistique de la composante du moment dipolaire d’une 
molécule selon l’axe Z (dans la direction du champ extérieur) s’exprime par la formule 


[mre exp | HR dQ 


| e یں‎ dQ 


où e = (E)! E et dQ représente langle solide élémentaire. 

Dans le cas où les molécules présentent une symétrie de révolution, z étant l’axe 
de révolution, u, = Ug = 0, ویر‎ = U; % = Xa = Ë j as = a; Qa = 0 pour AF ولا‎ 
Qu = Qa = — $ Qs3 = — Q (cf. Böttcher 1952, § 13), où Q équivaut à $ eq de Zwanzig 
(1956). 

En développant l'expression (21) en série selon les puissances de l'intensité du 
champ E,, et en se limitant aux termes contenant E,, on trouve expression suivante: 

u? E 16 Qcs E ate 
سو ان‎ VETO Te a) Boa tT- fat * 


)22( | 7 سیسات ++ 


<m -eò = (21) 


رو ژر- 1 ۰ ۰ج وکر- 1 
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et.en écrivant que 


3(e — 1) ot a+ | 8۱ : 2a) ) =z 
00 77 


Co S = 


on déduit pour la polarisation par unité de volume 


e—l N p? 
ERSE pya E 
2 47 3:7 (1 - fay)? ° 
pe ات‎ Le (23) 


a "Be + 2) at 090 ۶ہ ر-‎ 


4, Discussion des résultats 


Pour les particules n’ayant pas de moments quadrupolaires, on peut déduire de 

l'équation (23) à la suite des simplifications usuelles, la formule bien-connue de Onsager 
(1936). Une comparaison des valeurs de yw pour les liquides, obtenues avec la formule 
de Onsager, et de celles pour les mêmes substances en état gazeux a été donnée par 
Böttcher (1939 où 1952), et ensuite corrigée et élargie par Smyth (1955, tableaux 16.1 
et 16.2). La formule (23) pourrait bien servir à corriger la valeur de y pour les liquides, le 
moment quadrupolaire étant connu par une voie différente. Nous ne connaissons pas 
actuellement les valeurs des moments quadrupolaires des substances énumérées par . 
Smyth; d’autre part, pour quelques substances dipolaires peu nombreuses ayant des 
moments: quadrupolaires connus (Smith et Howard 1950, Gordy, Smith 1955 tableau 
31) l’auteur n’a pas disposé de toutes les données nécessaires pour les calculs analo- 
` gues à ceux de Böttcher (1939). 
On peut pourtant évaluer le deuxième terme “quadrupolaire“ de la formule (23) 
d’une manière approximative. Pour le moment quadrupolaire de quelques 10-26 u. e.s. 
et pour des valeurs ordinaires des autres paramètres, la valeur de ce terme peut atteindre 
d’un à vingt pourcent du premier terme “purement dipolaire“. Cette valeur dépend 
d’ailleurs très fortement de celle du “rayon de la molécule“ a 

Les termes quadrupolaires contenant (E,)? ont été تو‎ calculés. Toutefois, 
ceux-ci ne présentent pas d’importance pour la saturation diélectrique, étant de valeurs 
trés réduites. 

L’équation (23) a la forme P = Au? + Bu+ C; pour des particules n’ayant pas de 
moments quadrupolaires, celle-ci prend la forme P = Au?+ C. Pour les molécules 
présentant un moment quadrupolaire Q < 0, B < 0 aussi, et la valeur de مم‎ donnée 
par (23), la polarisation étant donnée par l’expérience, dépasse celle qui résulterait ` 
pour la même substance considérée comme dépourvue de quadrupôles (c’est-à-dire 
en utilisant la formule de اس پا‎ C’est le cas des dérivés halogénés des ny de 


1 Dans ce dernier ouvrage on trouvera une bibliographie du sujet allant jusqu'a 1954. 
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carbures, et les tableaux de Smyth (1955) montrent que la valeur du moment dipolaire 
permanent calculée avec (23) sera plus proche de celle pour l’état gazeux, car la valeur 


¢ c6 


gaz“ est plus grande que la valeur “Onsager“. 

On remarquera que ces derniers composés possèdent un moment Q < 0 en con- 
sidérant une molécule à deux ions comme representée par un modèle naif de deux 
sphères portant des charges électriques à signes opposés; si la plus grande sphère 
porte la charge électrique positive, le système a un moment quadrupolaire négatif, 
et inversement. Pour les hydrocarbures deux fois halogénés, évidemment Q < 0. 
Pour les autres composés de ces tableaux, par exemple les cétons, il faudrait procéder 
à une analyse plus approfondie du signe de leur moment quadrupolaire. 


` 


L’auteur tient à exprimer sa reconnaissance à M. le Professeur S. Szczeniowski 
pour de nombreuses discussions qui ont facilité la préparation du présent travail. 
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INVESTIGATIONS OF BREMSSTRAHLUNG OF ELECTRONS 
IN THE ENERGY INTERVAL 102 — 10% eV 


By J. BENISZ, Z. CHYLINSKI AND W. WOLTER 
Institute of Nuclear Research, Kraków 
(Received July 19, 1958) 


Four high energy (~1012 eV) electron-photon cascades have been investigated at 
the first stage of their development. The number and the energy of the electron pairs of 
the first generation, produced on the first radiation length were estimated and compared 
with the theoretical values calculated on the basis of the theory of Bethe and Heitler. 
On the other hand, the same was calculated according to the theories of Landau, Pome- 
ranchuk and Ter-Mikaelyan which take into account the influence of the medium on the 
bremsstrahlung of electrons of very high energy. The experimental results are in better 
agreement with the predictions of Landau, Pomeranchuk and Ter-Mikaelyan than with 
that of Bethe and Heitler. This fact confirms the results obtained earlier in our laboratory 
(Miesowicz et al. 1957). j 


Introduction 


Up to the present time there were very few experiments investigating the problem 
of the influence of a dense medium on the probability of emission of bremsstrahlung 
photons by electrons of very high energies, following the ideas suggested by Landau, 
Pomeranchuk and Ter-Mikaelyan. According to the theories developed by these 
authors we ought to expect that for very high energies (> 101? eV in nuclear emulsions) 
the spectrum of bremsstrahlung photons is strongly cut down for low energy photons 
as compared with the Bethe-Heitler spectrum. 

The existence of the effect predicted by the theories of Landau, Pomeranchuk 
and Ter-Mikaelyan (which will be further denoted by L.P.T.), was reported in a previous 
paper from our laboratory (Miesowicz et al. 1957). In that paper an electron-photon 
cascade generated by a photon of energy —101? eV has been analysed and a deficiency 
of low energy photons in comparison with cascade theory has been found. It has 
seemed that the observed lack of low energy photons was outside the limits of 
probable fluctuations. 

For studying the quantitative agreement of the observed effect with the theory, 
more cases of very high energy electron-photon cascades were needed, for increasing 


(143) 
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the statistical significance of the obtained results. However, electron-photon cascades 
of energies ~10!2 eV are very rare in photographic emulsions. 

In our laboratory in connection with the investigation of high energy nuclear 
interactions (jets), several liters of emulsion have been scanned for electromagnetic 
cascades. The aim of this paper is to present an analysis performed on four electro- 
magnetic cascades with primary energies 1012ھ‎ eV found in this scanning. The 
results of these measurements were compared with the theoretical predictions of 
Bethe and Heitler (B. H) and L.P.T. 

In the last few weeks after closing our work we received a preprint by Varfolomeev 
et all. in which the discrepancy between the observed development of an electron- 
-photon cascade and that predicted by cascade theory, has been also interpreted on 
the basis of the L.P.T. theory. 

Landau, Pomeranchuk (1953), Ter-Mikaelyan (1953, 1954) (Feinberg 1956) 
have called in question the validity of the B. H. formulas for very high electron energies 
and for photons for which W/E <1, where W denotes the energy of the emitted 
photon and E the energy of the radiating electron. In this case the B. H. formulae 
are distinctly different from that, given by Landau, Pomeranchuk and Ter-Mikaelyan. 
According to L.P.T. the reasons for this difference are as follows. If W/E is sufficiently 
small (W/E <1) the change of the momentum of an electron in consequence of 
a photon emission, e. g. in the Coulomb field of a third body, is so small that the 
uncertainty of the localization of the event, due to Heisenberg’s principle, is much 
greater than the distances between the scattering centres. If so, in Landau’s opinion, 
the effects due to the particular centres cannot be treated additively as assumed in 
the B—H calculations. The passage of an electron through the medium must be treated 
integrally. In the first estimation, Landau assumed the validity of the classical electro- 
dynamics in its semi-classical treatment which is justified by the high energy values 
of the radiating electron. In these calculations the process of emission is conditioned 
by the entire motion of the electron. The formulae obtained do not essentially differ 
from those, obtained later on the basis of quantum mechanics by Migdal (1954, 1955, 
1957). 

The following formulae are taken from the paper of Feinberg (1956). If the energy 
of the primary electron E is greater than ہیظ‎ where Eri پک سا‎ (=a f 
x me ~9 x 10° eV in nuclear emulsion E, = (4x M mc? = 21 MeV 
L is the radiation length in cms equal to 2.9 cm in nuclear emulsion, and n is the 
number of atoms in 1 cm? in a given medium, then the energy interval of the photons, 
between the energies 0 and E may be divided into three parts: 

(1) For photons of energies 


[ol 2 \‏ ہ١‏ مت 
WW, > 7 (422264) ous ) = 0.16 E% eV‏ 
m. cE?‏ 


1 We are much indebted to Professor Gurevich for sending us the preprint before publication. 
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the intensity of radiation is modified by the polarization of the medium. This effect 
has been evaluated by Ter-Mikaelyan 
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(2) For W, <W <F, 


a Ne g‏ او 
a + = 4.3 x 10-4 Et eV‏ ) ]ت۶۱ 


we have the modification of Landau and Pomeranchuk due to the effect of multiple 


scattering, 
es E. wm \” 
رو و ا‎ s 


(3) For energy values greater than W,, the B. H. formula remains in force: 


2 2 
dtan = = (| ر‎ 


3r \mc?) L 


For example, for the primary energy E = 1012 eV we obtain by means of these 
. formulae for nuclear emulsions the following values: W, = 1.6 x 107 eV and W, = 
= 4.3 x10! eV. From these values, we see that the modification which results from 
the polarization effect of Ter-Mikaelyan is negligible, thus the main modiflcation is 
due to the effect of multiple scattering evaluated by Landau and Pomeranchuk. 

In our case, we have calculated the average number of electron pairs of the first 
generation with energies greater than W which are created by the primary electron 
of initial energy E,(E, = [E(t)],-9) on the length t, where t is the distance from the 
beginning of the electron trajectory. The calculations were made, first on the basis 
of the B. H. formulae and then also by means of the L.P.T. formulae. 

Fig. 1 shows the integral energy spectrum of electron pairs of the first generation 
obtained from these two theories. 

We have taken into account the energy loss of the primary electron resulting from 
radiation. The energy losses resulting from ionization (about 2 X10? eV on one radia- 
tion length) are in our case negligible in comparison with the energy of about 101? eV. 

In our problem we must take into consideration the cross-sections of the following 
processes which are decisive in forming the energy spectrum of the electron pairs. 


(1) Pair production by photons 

(2) Compton effect 

(3) Photonuclear reactions 

(4) Production of tridents. 

Processes (1), (2) and (3) are competitive. One can show, and it results from 
` the geometry of the experiment, i. e. from the fact that the number of pairs is observed 
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on the first radiation length, that the reduction of the number of pairs caused by (2) 
and (3) is negligible. The number of pairs of energy W > 107 eV is diminished by 
about 2—3 %. The number of pairs of energy W > 108 eV rests practically unchanged, 
because in this energy range the cross sections of processes (2) and (3) are very small 
in comparison with the cross section of (1). 

Increasing of the mean free path for pair production with decreasing energy of 
the photon (10? eV) works in the same direction. Although both causes act in the 
same direction as the modification of L.P.T, nevertheles their effectss are much smaller. 
We have estimated that their whole cumulative influence amounts only to same per 
cent and so, it is negligible in comparison with the L.P.T. modification. 

Thus we see that the number of electron pairs of the first generation is a sensitive 
detector of the character of the bremsstrahlung spectrum and quite insensitive to com- 
petitive effects, in comparison with the process of pair production. 

The number of pairs is of course increased by the trident production, however, 
this does not privilege the production of pairs of small energies. Accepting that for 
E, = 5X10" eV, the mean free path for trident production is about 3.5 rad. lengths, 
we obtain a contribution of about 6% of the whole number of pairs produced by the 
first process. Finally we see that the effects (2), (3), and (4) are of no importance in 
our problem. 

On Fig. 1, representing the results of both theories (B. H. and L.P.T.), it is 
interesting to remark that in the case of the B H. energy spectrum we have a depen- 
dence of the number of pairs on E,/W only, while in the analogous curves of L.P.T. 
there is a dependence on both E, and W. 

Although there is a possibility of large errors in estimating Ey = the primary 
energy of the electron, it is important to stress that even a change of E, of one order 
of magnitude, does not change the character of the L.P.T. spectrum. From Fig. 1 
we see also that the L.P.T. spectrum curves have a plateau beginning from a certain e, 
which feature is not shown by the B.H. spectrum curve. This is the main consequence 
of the L.P.T. theory which just expresses the lack of small energy pairs. 


Experimental Method 


In order to compare the experimental data of the energy spectrum of brems- 
strahlung photons radiated by an electron of high energy with the theoretical predic- 
tions of B.H. and L.P.T. we have adopted the following procedure: 

(1) We have evaluated the energy of the primary electron pair from the develop- 
ment of the electron-photon cascade. 

(2) We have evaluated the energy of all electron pairs, produced by the conversion 
processes of bremsstrahlung photons, emitted by electrons of the primary pair. It 
is obvious that the distance from the point of origin of the primary pair up to the 


end point where the scanning for electron pairs is stopped ought to be as long as 
possible. 
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The restriction to electron pairs coming from the first generation only, was imposed 
by the conditions of our measurements. 

Three of the investigated cascades originate as a result of decays of x° generated 
in the same nuclear interaction, therefore the mutual radial distances between the sepa- 


rate cascades should be small. This fact does not allow to conelate with the particular not 
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Fig. 1. Integral energy spectrum of electron pairs of the first generation. 
N, — number of pairs of energy greater than the given value £. 
Full curves — Bethé-Heitler curves for three different depths (1.0; 1.5; 2.0 rad length). 
Dotted and dotted-dashed curves — Landau — Pomeranchuk and Ter-Mikaelyan curves for the depths 
1.0, 1.5, 2.0 rad. length and the primary electron energies Ej = 3 X 1011 eV; 1 x 1072 eV and 3X10! eV. 


cascade any electron pair generated in a greater distance from the tracks of the primary 
electron pair. On the other hand, it is known that all the electron pairs of the first 
generation, in consequence of the great energy of the primary electrons, must originate 
in close proximity of the parent track as apparent tridents. Since the losing of such 
pairs is rather improbable, our procedure releases us from the scanning bias. 
The experimental material presented in Table 1 contains four electron-photon 
cascades. Under the assumption of equipartition, we attributed to the six primary 
electrons of the three electron-photon cascades equal energies of 6.0 x 101" eV. The two 
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Table 1 


Bees From angle of 
£ z تج‎ From scattering measurements divergence 
Pair TELA RES Remarks 
$ 5 ee = E+ وط‎ 
وق‎ E, (eV) E, (eV) (ev) W (eV) 


Cascade A 

Primary pair *£,+ 
+E,=1.2X102 eV 
7.0 X108 1.2 x 109 2.0 x 109 First generation 


1.3 x 108 6.6 x 108 7.9 x 108 Order of generation 
impossible to esta- 


blish 
>2.4x 09 >2.5 x 109 >4.9 x 109 — First generation 
= = = +10 x10? First generation 


(3.0*2-8) x 108] (3.579) x 108] 6) x 108 — Order of generation 


impossible to esta- 


blish 
>4.6 x 109 First generation 


First generation‏ ۶ ھپ تی 


Primary pair *E, + 
+ E,~1.2 1012 eV. 


(7.8+15:7)x 108|- 0ن‎ SE First generation 
>3.4x 109 > 3,8 x 109 57.2 X10? First generation 
>3.5x 109 > 3.6 × 10° First generation 


— First generation 


First generation 


1 Primary pair *E,+ 
+ نہ و‎ 1.2×1018 eV 
2 LHC و‎ (1.3*2-6) x 109] (1.4#2-7) x 109] 25x10) First generation 


3 >3.2 x 10° >4.8 x 109 >8.0x 10? >1 101° First generation 


72 


* Eı +E, = 2E, according to the assumption of equipartition of energy. 
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Ed. Table 1. 


ESET From angle of 
as 7 E From scattering measurements di 
Don lvergence 
Pair |4" 3 _ À 
Fas سی تر و‎ Remarks 
اڈ وق‎ E, (eV) E, (eV) meee! 
AEDES (eV) 


0.57 | (1.375'3) x 109 (3.03 :8) x 108| (4.3+4:1) x 108 First generation 


0.89 (2.273?) x 108 (2756) x 108 (4.917) x 108 Order of generation 
impossible to esta- 


blish 


Cascade D 

= — — Primary pair *E,+ 

+ E,~7.0X10!! eV 
2 0.32 2.0 x 108 1.0 x 108 3.0 x 108 — First generation 
3 0.57 1.3 x 108 1.0 7 1.4.x 108 — First generation 
4 0.60 6.0 x 108 = ۶ 08 — First generation 
5 0.82 1.3 x 108 8.0 x 107 2.1 x 108 — First generation 
6 0.96 >1.0x10? 7.5 × 108 > 18x10? — First generation 


remaining primary electrons of the fourth cascade have both the same energy 3.5 x 10" 
eV. In order to make possible a common treatment of all cascades, the experimental 
material was compared with the curves of B. H. and L.P.T. calculated for the primary 
electron energy E, equal to the mean energy value of all eight primary electrons. 

Three cascades (A, B, C) (see Table 1) have been generated in the same nucleon- 
-nucleon interaction of the type 0+14g (Ciok et al. 1957). The energy of the primary 
particle was found to be 3.3 x 1014 eV/nucleon as obtained from the angular distribution 
of shower particles. The cascades, mentioned above, are contained within the narrow 
cone with the angle of the order of 10.4 radians with respect to the direction of the 
primary particle. The distances from the point of interaction to the point of origin of each, 
of the four primary pairs are respectively: 0.97, 1.17, 1.28, 1.58 radiation lengths. 
Up to the depth of about 2.2 radiation lengths from the point of interaction no more 
primary electron pairs were found. Since at greater depths than about one radiatior 
length from the point of origin of each primary pair, electron tracks belonging to 
different cascades are crossing each other, there was no possibility of evaluating the 
energy of each cascade separately. For this reason we evaluated the mean energy of 
the cascade, under the assumption that there is energy equipartition as well between 
the primary photons as between the electrons of the primary pair. The mean energy 
of the cascade has been evaluated from the longitudinal development of all four 
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cascades assumed to be one cascade initiated by four primary electron pairs of equal 
energy, created at different depths. A correction for lateral distribution of the cascade 
has been taken into account. 

A target diagram at the depth of 2.76 radiation lengths from the point of inter- 
action has been made. We have received the following data for the number of electron 


tracks: 


59 electrons of energy greater than 1X10°eV within a radius of 80 y and 
97 electrons of energy greater than 410% eV within a radius of 280 مم‎ 


The energy values have been determined by means of scattering measurements. After 
taking into account the correction for lateral distribution according to Pinkau (1956) 
we received the following data: 


(a) 70 electrons of energy greater than 1 X10° eV 
(b) 101 electrons of energy greater than 410% eV 


From (a) we have obtained for the mean energy of each of the four primary photons 
Eon = (1.1794) x 1012 eV and from (b) E,, = (1.3493) x 10!%eV. In these calcu- 
lations we used the cascade tables of Janossy (1951) (electrons produced by electron) 
as well as Janossy’s (1952) standard deviation. 


In the sequel we have accepted the value 1.2 x101? eV for the mean energy of 
each of the four primary photons. The mean energy of the 2° meson evaluated in 
this manner (assuming equipartition of the energy of 2.41012 eV between the two 
photons) is in good agreement with the energy values of the two charged mesons, 
contained in the same angle as the cascades. Their energies, obtained from an analysis 
of secondary interactions are equal to 1 X1012 and 2 × 1012 eV respectively. 


Now we have tried in turn to evaluate the energies of the electron pairs of the 
first generation for each cascade separately. The correlation of a given pair to the 
first generation pairs was done by using the criterion of energy and the criterion of 
geometry. Such a procedure was successful in all cases except in the casé of three 
electron pairs where there was no possibility of establishing the order of generation, 
i. e., in which succesive generationthe pair was produced. Asa matter of fact such a pro- 
cedure could be adopted in only three cascades out of the four registered. In the case 
of one cascade neither the number of the first generation pairs nor their energies 
. could be evaluated. 


The energy measurements have been performed by scattering, especially the 
differential one. In cases in which scattering measurements were not possible, we 
have based on the data derived from measurements of the angle of divergence of the 
two tracks of a pair. 


In Table 1 are given the results of our measurements performed on four cascades. 
The data concerning cascade D initiated by a single photon of energy Ean = 
= (7.0*55) x 1011 eV were taken from the paper of Miesowicz et al. (1957). 


The energy spectrum of electron pairs of the first generation created on the 
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first radiation length has been done on the basis of the data given in Table 1 (Fig. 2). 
The histograms represent the experimental spectrum taking into account the uncer- 
tainty of the order of generation of the three electron pairs (histogram a and b). For 
energies greater than about 10° eV (see the point with an arrow on Fig. 2), the shapes 
of the histograms are unknown because of the impossibility of evaluating the upper 
limit of such high energies. Curves I and II represent the B.H and L.P.T. spectrum 
respectively, for primary electron energy equal to 5 × 1011 eV at the depth of one 


radiation length. This energy value corresponds to the mean energy of eight primary 


electrons under consideration. The dashed curves represent Poisson’s deviation from 


Np 
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Fig. 2. Integral energy spectrum of electron pairs of the first generation 0 on the first radiation length. 
Np — number of pairs of energy greater than the given value W. 
a, b — experimental histograms. | 
I, II — Bethe-Heitler and Landau — Pomeranchuk —- Ter-Mikaelyan curves for the primary electron 
energy 5X10! eV respectively. 
Curves I and II are given with their standard deviations. 
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the theoretical curves of B.H and L.P.T. It is worth while to remark that taking into 
account the uncertainties in the energy estimation of the electron pairs (about 100%) 
does not influence to a high degree the shapes of the histograms. 


Conclusions 


1. The experimental energy spectrum of the electron pairs of the first generation, 
produced on the first radiation length, shows a statistically significant deviation from 
the Bethe-Heitler energy-spectrum curve. 

2. There is rather a good agreement between the experimental results and the 
curve, which represents the energy spectrum of Landau, Pomeranchuk and Ter- 
-Mikaelyan. 

3. The method of investigation of the electron pairs of the first generation only, 
is in the opinion of the authors a sensitive tool in detecting the difference between 
the energy spectrum of Bethe-Heitler and Landau, Pomeranchuk, Ter-Mikaelyan, 
since in the cascade development there is a degradation of energy of the emitting 
electron. 

We wish to express our gratitude to Professor M. Miesowicz for suggesting these 
investigations, his permanent interest in our work, and many valuable ideas and hints 
during stimulating discussions. 
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The statistical theory of order-disorder transformations for n-component alloys with 

` arbitratry composition and crystalline structure is given. The parameters of long-range 

and short-range order are defined and equations which hold for them are found. The 

configurational free energy in zeroth and first approximations is calculated on the basis 

of generalized Bragg-Williams and quasi-chemical methods. By taking in to account 

the interaction of atoms from more distant shells than the first one it is passible to compute 
more precisely the order-disorder transition temperature for an arbitrary alloy. 


1. Introduction 


Several theoretical methods of investigations of order-disorder transitions in 
binary alloys are being elaborated (Guggenheim 1952, p. 101). Though these theories 
are in good agreement with experiment, they do not allow to treat uniformly order- 
-disorder transitions in alloys of arbitary composition and of various crystalline structure. 
Moreover the theory of order in more than two component alloys was not elaborated 
till now. It is therefore the aim of the present paper to give a general theory of order- 
-disorder transformations for n-component alloys with arbitrars composition and 
structure n-component alloys, the superlattices of these alloys consisting each of l 
sublattices. 


2. Long-range order parameters 


In the general case of an n-component alloy, its superlattice may be examined as, 
‘a superposition of / mutually interpenetrating sublattices. Each of these sublattices 
has M, lattice-sites (u = 1,2, ..., I). We denote by N the whole number of atoms 
in the alloy, by N; the number of atoms of kind i (= 1, 2, ..., n), by a;, the number 
of atoms of kind i in the y-th sublattice, and by p,, the a priori probability that 
an atom of kind i is in any site of the y-th sublattice. 
(153) 
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For p;, we have the relation 
ais (1) 
Piu M, : 
We introduce m linearly independent long-range order parameters s“ (k = 
= l, 2, ..., m) by defining them as linear functions of p;, and requiring that 


__ ÍL, for the state of perfect order, 
°k = ١0, for the state of perfect disorder. 


Since for several alloys we know their superlattices, we can determine pf,, where pj, 
is the p,, for the state of perfect order. Defining the perfect disorder as a configuration 


of atoms for which p;, = ¢ = ae we can define s, generally and uniquely by 
writing the equation of the straight line passing through the points (1, p§,) and (0, c;); 
thus 


Se = Sp = PE, = retê (2) 
Pipe ب‎ 
The number m appearing above determines the number of long-range order 


parameters which must be introduced for describing the configuration of atoms 
in a n-component alloy (Wojciechowski 1958a). It is easy to show that 


m=nl—n—1+1. (3) 


3. The zeroth approximation 


In the zeroth approximation we may neglect the correlation of atoms in the 
lattice-sites and treat a,, as independent quantities. Denoting by E(s,) the average 
configurational energy of the alloy, we may write for the configurational free energy 
F(s,) the following expression 


F(s) = E(s) — kT In W (5p); | (4) 


where W (s,) denotes the number of ways of arranging a,, atoms over all lattice-sites 


7 )+۶ 080 9 ۷ ئ) 
1 
W (s,) = [I W, (sp) (5)‏ 
u=1‏ 


where 


M,! 


n 


Wp (sk) = 


di ! (6) 
A K 


t 


In order to get an expression for the average configurational energy for ae 
n-component alloy of arbitrary crystalline structure and composition, we shall assumn 
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that this energy depends upon the ordering of atoms in the alloy, and that it is a super- 

position of interactional energies of these atoms. Thus, denoting by V,.(0) the inter- 
۰ . . J 

actional energy of an atom of kind k with an atom of kind J, their mutual distance 

being ,و‎ and by q,,(@) the average number of pairs of atoms E —j with mutual distance ہو‎ 

we may express the total configurational energy, as follows 


qe (0) Va; (0). (7)‏ 2 وو اش 
It is obvious that ils‏ 
Ge (0). (8)‏ = (وایرو Vo) = Vo),‏ 


Since each lattice-site of the u-th sublattice has Zu» (0) neighbouring lattice-sites 
in the »-th sublattice, their mutual distance being ọ, we get the average number of 
pairs ,;(0) by summing up the products z,, (0) a, Pj over all lattice-sites and dividing 
the result by 2. Therefore 


I 
dej (o) = 3 Dj Zur (o) Aku ٢۰ (9) 
Hv 


Using Eq. (9) and taking into consideration that a;, = f;,(s,) and pj, = Eiu (Sh) 
we may rewrite Eq. (7) in the following form 


2001 
E (sy) =— 3 DD y Z uv (o) V; (o) Qiu (s4) Pjv ٭(یرد)‎ )10( 
e ij m» 
In the papers dealing with binary alloys “the nearest-neighbours“ assumption 
(o = 1) has been made, thus from (7), we have 


E=- } X yO yA. (11) 


It follows therefore from Eq. (10) that the above assumption is not essential in our 
theory. Furthermore the approximation of Eq. (10) consists in neglecting the corre- 
lation. To take into account this correlation, we introduce the quantities Q% (o) denoting 
the numbers of pairs consisting of an atom of kind A in the w-th sublattice and an atom 
of kind j in the »-th sublattice, their mutual distance being .و‎ In the next part of 
this paper we express by Q47(0) the short-range order parameters and the configura- 
tional energy of the alloy. 

Substituting Eqs. (5) and (6) into Eq. (4) and using Stirling’s formula, we find 


l n 
F(s,) = E(s;) =F kT D, > all a, — 1) + const. (12) 


u=1 i=1 


` The dependences s, = f,(T,V) may be found from the equilibrium conditions 


1 n 
OF (sk) _ OE(s%) ota 1 
SS er DD In ay, = 0. (13) 


= 


156 Kazimierz F. Wojciechowski 


4. Short-range order parameters 


Summing Qs; (g) over all k,j and u,v we get the equality 


lon ہر ھا‎ ln n:n 
D 3J Qs (e) +B {DJ 3) 0 (e) + DUH ()م)‎ = yo (14 
v=1k=1 u*v k=1 pv» ky k=1j=1 


in the above equation the factor 1/2 appears since Qe; (e) = Qj (0). On a bassis of eq. (9) 
we find 


n n PAT 
> 2 dki (o) = i> 2 Dies (0) Aku Piv: ۱ )15( 


Omitting in the eq. (14) the summation over y and k, and taking under consideration 
(15) for case when ¢ = 1 and Qg (1) = 0, we may rewrite equation (14) in the 
form 


1 n 1 
2 {Où OR 2 QO (D) = au È w (1). (16) 


If we restrict ourselves to the above case than (16) represents n+l constraining 
relations between Q4;(1) and a,,. An analogous dependence can be obtained in the 
general case using (14). In both cases the number of pairs Q% (o) is greater than the 
number of relations among them. Thus we have to introduce additionaly, by definitions, 
as many short-range order paramaters as are necessary to express all Qj; (o) by them. 
If we denot by Q; (9) the independent Qf; (g)appearing in the equation (16), or in the 
corresponding general equation, then we define i short-range order parameters رہ‎ (0) 
as follows 


Q; (e)/Q(@) - 2c J Fes 


QIQ) = 2ee cj 00 


| = وہہ = وا cu‏ 


where Q? (0) denotes Q; (o) in the state of the complete order and Q(o) denotes 
the whole number of pairs for given (g). It is easy to compute Q? (0) for each 
alloy. 


5. The first approximation 
In the first approximation we take into account the correlation of pairs of atoms 


generalising the quasi-chemical method. Now the configurational free energy can be 
expressed as follows 


F (sp, وہ‎ = U (4,5) - دا7‎ W (54,0), ` | (18) 
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where U (s,, 0;) is connected with the energy V, of state t by the following relation: 
D e-ValkT - y (sy, 8) e— U (sp oj)IRT, 
lEn 
a 22 Qk; (0) Vay (0). (19) 
ts» Rs j 
ںہ‎ is the most probable value of o; and is given by relations 
9 x o 
a {In W (s,,0;) — E (sg, o;)/kT} = 0. l (20) 
Now E (sg, o;) = 2 (s,,6;) and U(s,, ہہ‎ can be calculated from the relation 
9 
1 
(7) 
Therefore we see that now the problem now reduces to the calculation of W (sg, 0;). 
To compute W (sg, 0;) we shall apply the quasi-chemical method assuming that 


W (sg, o;) is proportional to the total number of ways one can divide Q (e) pairs into 
groups of pairs Q%; (0). Thus 


E (Sk, Gi) = 


[U (sr ۰ء‎ (21) 


ITQ (0)! 
Kise رہ‎ KG a کے‎ ` (22) 

IITTI پان‎ (o)! 

e Hvk,j" 
The proportionality factor K(s,) appears in the eq. (22) since J7Q (0)!/IT IT IT Q% (o)! 
contains many unreal situations which occur owing to the existence of correlation 
between the pairs Qf; (0). In order to determine K (s) we shall consider of which 
maximize the function W (sg, 0;). As it is known from statistical mechanics, the maxi- 
mum value W (s,, of) approximates very closely the sum of W (sg, 0;) over all o; 
for a great many assemblies. Therefore 


W (sp, of) = 2 W (sg, 0) = W (s), (23) 


where W(s,) is given by eq. (5). 
Putting in eq. (22) o; = of and applying eq. (23) we find 
i lon 
IT IT IT Q% (e, اہ‎ 
K = W @ MY Ro] 
(sai HOE TT IT 


where Qi (0, of) denotes Q¢; (¢) in, the state of complete disorder, and 


| 
Qf (e, of) = 5 Zur (2) aru Pjr (232) 
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Finally we get 


lon 
IT TI Q% (e, اہ‎ 
W (sk, o) = W (s) I F. (24) 
3۳/7100 
uv kj 
Since we can express all Q% (e) by o; and all o; by s, (by (20)) we obtain finally 
the configurational free energy sought for. We can represent this energy for an arbitrary 
n-component alloy as a function of long-range order parameters only. 


6. The influenece of interaction of farer neighbours 


No theory of order in binary alloy elaborated thus far has been able to take into 
account the influence of the interaction of farer atoms farther away, than the nearest 
ones. The present theory allows permits the treatment cf n-component alloys of 
arbitrary with respect to crystallic structure and composition. It allow also makes it 
possible to take into account the influence of the interaction of atoms from shells 
more distant than the first one. 


Fig. 1 Fig. 2 


Fig. 1. Simple cubic structure a denotes lattice constant. 


_ Fig. 2. Body centered cubic. structure. © denotes lattice sites of the sublattice 1; @ denotes lattice sites 
of the sublattice 2. 


The quantities characterising the given superlattice are: the distances o of atoms 
in the lattice of the alloy from a definite atom in this alloy, and the coordinational 
numbers z,, (0). In the table I are given the succesive distance o = 192-5 Sy for 
a simple (Fig. 1), body-centred (Fig. 2) and face-centred (Fig. 3) cubic structure. These 
distances are denoted on the figures and in the table by 0(1), @ (2), ..., 0(5). 

From the table. it readily follows that for a simple cubic structure e(1), e(2) 
and 9(3) are magnitudes of the same order. But for the body-centred and face-centred 
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Table I 
0 Simple cubic Body c. cubic Face c. cubic 
av3 2 
e(1) | a a F 0.86602a K 0.70710a 
o(2) aV 2. 1.414218 a a a a 
0(3) | ۷3 1.732054 a2 1.41421a aV2 1.414160 
۵ہ‎ | 2a 2a av3 1.73205a av3 1,73205a 
e(5) کے‎ — 2a 2a 2a 2a 


cubic structure the distances p(2), o(3), o(4) being magnitudes of the same order, 
differ from o(1) by one degree of order. 

The interactional energies پر7‎ (0) depend strongly on ọ thus if we wisch to obtain 
better results than that obtained from the “nearest neighbours assumption“ one has 


Fig. 3. Face centred cubic structure. © denotes lattice sites of the sublattice 1; @ denotes lattice sites 
of the sublattice 2. ` 


to take into account the influence of atoms from second, third and fourth shells has 
to be taken in to account. 

This allows, for instance, a more precise determination of the order-disorder 
transition temperature, which for an arbitrary alloy, as pointed out by Swierczkowski 


and Wojciechowski (1958), depends only on بر‎ (o). 


7. Conclusions 


The above theary allows us to consider of long and short-range order in n-compo- 
nent alloys of arbitrary crystallic structure and composition. The zero approximation 
of this theory allows, by regarding the influence of atoms from an arbitrary shell, 
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a more precise determination of the order-disorder transition temperature specific, 
heat, etc. The only parameters of the theory are the interactional energies V}, (0) 
which are to be determined by experiment. 

The first approximation allows us to take into account the correlation in arbitrary 
alloys, but as it has been mentioned in the papers, dealing with the ternary alloys 
(Wojciechowski 1958a, b) in alloys of more than 2-components the first approximation 
leads to complicated mathematical problems whose exect solution seems to be very 
difficult. Therefore we may suppose that for practical purposes the zeroth approxi- 
mation is the only suitable one. 

Examples of applications of this theory to ternary alloys are given by the author 
in (1958a, b). 

The author wishes to express his gratitude to Professor J. Lopuszañski for helpful 
advice and discussions. 
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It is shown in this work that an LC resonant circuit in the anode circuit of a vacuum 
tube controlled by a piezoelectric quartz crystal cannot be used unrestrictedly for the 
measurement of capacity and dielectric constants. This method is limited by the dielectric 
losses and conductivity of the condenser under measurement. 


It is shown in this work that an LC resonant circuit in the anode circuit of a vacuum 
tube controlled by a piezoelectric quartz crystal cannot be used unrestrictedly. for 
the measurement of capacity and dielectric constants. This method is limited by the 
dielectric losses and conductivity of the condenser under measurement. 

One of the most accurate methods of measuring electric capacity is the method 
based on the observation of the change in the anode current of a vacuum tube oscillator 
controlled by a piezoelectric quartz crystal as the capacity is changed in this circuit 
(Maier 1954). The firm Kip & Zonen have even put on the market factory equipment 
based on this principle. The principle of this method is shown in Fig. 1. A circuit 
consisting of a capacity C and an inductance L is in the anode circuit of the crystal 
oscillator while the quartz crystal is in the grid circuit. Oscillations corresponding 
to the resonant frequency of the crystal are produced in the anode circuit. Their 
amplitude and anode current depend on the resonant frequency of the LC circuit. 
If the condenser C is variable and the LC circuit has a resonant frequency greater 
than the resonant frequency of the crystal, then, as the capacity C is increased, the 
anode current drops to a minimum, which occurs when the LC circuit is accurately tuned 
to the frequency of the crystal. Further increase in the capacity produces a rapid 


(161) 
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increase in the current, and finally, at some value of capacity the oscillations cease. 
The variation of the anode current with the capacity in the LC circuit is shown in Fig. 2. 
In order to measure an unknown capacity X we first set a capacity standard C at some 
point of the curve of the function J = J(C) and then connect in parallel an unknown 
capacity X. The oscillations then cease. We change the capacity C so that the current 
detector again indicates the same value of current as before the unknown capacity was 
connected. The decrease in the capacity C, AC gives the value of the unknown capacity. 
The accuracy of the measurement and the sensitivity of the method greatly depends 


SSS 


100 110 120 130 140 150 


Fig. 2 


on which point on the curve of the function I = 1(C) the reference point was chosen. 
It is most convenient to choose it on the steep side of the curve in the neighbourhood 
of the point where the oscillations cease. In Fig. 2 this is the point denoted by the 
letter P. In addition, the precision of construction of the standard condenser C, the 
sensitivity of the current detector, and-the stability of the anode voltage has a large 
influence on the sensitivity of the method. Under suitable conditions this method 
permits the detection of changes of capacity of the order of hundreds and even 
thousands of pF. If, however, the condenser being measured has dielectric losses, 
then after it is connected, not only the capacity of the circuit changes but also the 
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attenuation. There arises the question of whether the current corresponds to the same 
circuit capacity, that is, whether the method gives correct results in the case of 
a condenser with dielectric losses and what is the limiting value of losses which still give 
satisfactory results. In order to settle this question by an experimental method, the 


48 ۸ R 


و۸ 150 


1۸5 


Fig. 4 


Fig. 5 


circuit shown in Fig. 3 was built and modified so that the BRED detector was under 
low voltage and one plate of the condenser was earthed. In this way the variation of 
the capacity between the condenser and earth was avoided. A plot or the anode current 
vs the capacity of the oscillating LC circuit was made. Next, radio resistances were 
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connected in parallel with the condenser C and for each of them the above-mentioned 
curves were plotted. The capacity of the resistors were measured separately by a diffe- 
rent resonance method. The value was on the average 0.8 pF. Their inductance was 
negligibly small in comparison with the inductance of the LC circuit. In this way 
a plot of J = I(C) was obtained for resistors of various values over the range 20 kQ 
to 10 MQ. Fig. 4 shows the plots for resistances of 48 kQ, 150 k, and 1 MQ and 
a plot made without resistors. The curve shown in Fig. 5 has been prepared from 
anumber of measurements for different resistors. This curve expresses the dependence 
of the apparent capacity of the LC circuit on the resistance connected in parallel 
with the condenser C for a fixed value of anode current in the neighbourhood of the 
point at which the oscillations cease. It is seen from the curve that beginning with 
a resistance of about 400 kQ the effect of the dielectric losses and conductivity of the 
condenser on the accuracy of the measurement begins to decrease in the direction of 
decreasing conductivity. Analysis showed that for a resistance of 400 kQ and capacity 


of about 0.8 pF the dielectric loss is of the order of 5 x 107% at a frequency of 
200 KC/sec. The error in the measurement of the capacity was in this case about 1%. 
For larger capacities the loss would obviously be smaller. 

From this investigation it follows that the method described here cannot be used 
for measurements of the dielectric constant without restriction, in particular, it cannot 
be used for investigating substances with large conductivity. The condensor with the 
substance under investigation should not be a conductor with an ohmic resistance 
smaller than a certain limiting value, below which the systematic error of measurement 
sharply increases. The same applies to the dielectric losses. The standard condenser 
with the substance under investigation cannot have large dielectric losses. In using 
equipment of the above-described type for measurements of the dielectric constant 
one should first determine the limits of its applicability in an appropriate way. These 
limits are different for circuits operating with different crystals. 

The author wishes to express his gratitude to Professor M. Jezewski for suggesting 
the topic and for his aid in the course of this work. 
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